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Abstract 

We review our state of knowledge of coronal element abundance anomalies in the Sun 
and stars. We concentrate on the first ionization potential (FIP) effect observed in the solar 
corona and slow-speed wind, and in the coronae of solar-like dwarf stars, and the “inverse 
FIP” effect seen in the corona of stars of later spectral type; specifically M dwarfs. These 
effects relate to the enhancement or depletion, respectively, in coronal abundance with respect 
to photospheric values of elements with FIP below about 10 eV. They are interpreted in terms 
of the ponderomotive force due to the propagation and/or reflection of magnetohydrodynamic 
waves in the chromosphere. This acts on chromospheric ions, but not neutrals, and so can 
lead to ion-neutral fractionation. 

A detailed description of the model applied to closed magnetic loops, and to open field 
regions is given, accounting for the observed difference in solar FIP fractionation between the 
slow and fast wind. It is shown that such a model can also account for the observed depletion 
of helium in the solar wind. The helium depletion is sensitive to the chromospheric altitude 
where ion-neutral separation occurs, and the behavior of the helium abundance in the closed 
magnetic loop strongly suggests that the waves have a coronal origin. This, and other similar 
inferences may be expected to have a strong bearing on theories of solar coronal heating. 

Chromospheric waves originating from below as acoustic waves mode convert, mainly to 
fast mode waves, can also give rise to ion-neutral separation. Depending on the geometry of the 
magnetic field, this can result in FIP or Inverse FIP effects. We argue that such configurations 
are more likely to occur in later-type stars (known to have stronger field in any case), and 
that this explains the occurrence of the Inverse FIP effect in M dwarfs. We conclude with a 
discussion of possible directions for future work. 
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1 Introduction 


Working during the early years of solar UV and X-ray spectroscopy, Pottasch ( 19631) found evi¬ 
dence for significantly higher abundances of Mg, Si, and Fe in the low solar corona than in the 
photosphere, and concluded, somewhat reluctantly that “the chemical composition in the solar 
atmosphere differs from the photosphere to the corona”. These elements, Mg, Si, and Fe, are all 
elements with first ionization potential (FIP) less than 10 eV, now known to be routinely enhanced 
in abundance in the corona with respect to photospheric values, a phenomenon that has become 
known as the “FIP effect”. High FIP elements such as O, Ne, and He, have much smaller abundance 
enhancements, or even abundance depletions in the corona. Although the possibility of elemental 
fractionation between the solar photosphere and corona onl y re ally began to be taken seriously in 
the mi d 198 0s, with the publication of influential reviews bv lMeveil ( 1985al lblf. recognizing the work 
of lPottaschi (119631 ) almost fifty years ago makes the problem of understanding the FIP effect nearly 
as old as that of coronal heating. In fact, modern models of the effect to be discussed in detail 
below, in which the fractionation is driven by the ponderomotive force of Alfven waves, make an 
intimate connection between the abundance anomaly and coronal heating mechanisms, such that 
the FIP effect may yield several important insights into the nature of the latter. This is a theme 
we will develop throughout this review. 

Thesolar FIP effect manifests itself in several modes of observation. Spectrosco py, as pioneered 
by Pott aschl ( 1963 ) and extensively reviewed at intervals ov e r the last 20 years (e.g.. lFeldman . 1992 : 


Feldman and Lamina . 2000l : Feldman and Widine . 2003 : Saba . 1995t h reveals the composition 


primarily of the “closed loop” corona. Meveit ( 1985a 3) also considered elemental abundances 


measured in situ in the solar wind, and in solar energetic particle events, which flow out along 
open m agnetic field lines. O bservations with Ulysses, the first mission to fly ov er the solar polar 
regio ns ( Wenzel^^ef^iL , 1992) revealed FIP fractionation varying with wind speed ( Zurbuchen et all . 


19991 : von Stei ger et al. . 120001 1. Slow speed solar wind had abundances resembling those in the closed 


loop solar corona, whereas high speed wind from polar coronal holes had a much lower level of FIP 
fractionation. With the advent of the Solar and Heliospheric Observatory (SOHO), a more coherent 
observational picture began to emerge. Coronal holes, the source of the relatively unfractionated 
fast solar wind, were themselves shown to have similar abundances to the wind emanating from 
them. The lower latitude closed fiel d corona was also s hown to have FIP fractionated plasma, 
similar to the slow speed solar wind. Schmelz et al\ ( 2012 1 give a modern view of coronal element 
abundances derived from these various solar physics sources, coronal spectroscopy, solar wind and 
solar energetic particles. 

Astrophysical EUV and X-ray spectroscopy, made possible by the 1990’s launches of the 
Extreme Ultraviolet Explorer (EUVE), the Advanced Satellite for Cosmology and Astrophysics 
(ASCA), Chandra and XMM-Newton, allowed coronal abundances in stars to be measured for 
the first time. Again, nearly two decades after the first such observations, various observations 
seem to be falling into place. In this review, we will attempt to synthesize these two strands of 
observations, solar and stellar, into one complete picture of coronal element abundance anomalies. 

Following this introduction, in Section [2] we briefly review recent developments in solar pho¬ 
tospheric abundances, where improved spectroscopic data and the application of 3D radiation 
transfer calculations have ushered in a revised solar composition. Section [3] describes the various 
facets of the solar FIP effect, while Section [4] surveys stellar FIP and Inverse FIP effects. Section [5] 
describes early attempts to model the solar FIP effect. Section [6] lays out the model advocated 
here, where the ponderomotive force due to Alfven waves propagating through, or reflecting from, 
the chromosphere accelerates ions up or down, while leaving neutral atoms unaffected. Section 0 
describes the results of such a model, and discusses its important implications. Section [8] concludes 
with suggestion for future research directions. 


2 









































































2 Solar Photospheric Composition 


2.1 Review 

Any work on solar coronal abundance anomalies must begin with reviewing the photospheric 
composition. Despite its long history, the composition of the photosphere, also taken as a proxy 
for “cosmic abundances”, has undergone significant revisions in recen t years. We take as our 
starting point the composition review of lAnders and Grevessel ( 1989 1, this being the standard 
solar composition in use for the early studies of the FIP effect, and for many years the default 
abundance set impl emented by s pectra l fitting software in use in X-ray astronomy. This was 
updated bv iGrevesse and Sauvall ( 19981) . who revised downwards by a small amounts N and O 
(0.13 and 0.10 dex, respectively) and also Fe which moved from 7.67 (on a logarithmic scale where 
the abundance of H is 12) to 7.50, in agreement with the meteoritic value. This l ast modificat ion 
stem med from impro vemen ts in atomic data used to analyze solar spectrum (e.g., Holweger_efjrl 


199ll: Biemont et<jd\. 1 199 lh. and is largely supported by more modern analyzes ( Asplund et al. 
2000a b; Asplund , 2000 ; Bellot Rubio and Borrerol . 2002 ). 


The next major revision came to the photospheric abundance of O by lAllende Prieto et al\ 
( 20011 ). who applied a three-dimensional time-dependent hydrodynamical model solar atmosphere 
to the observed 6300 A line. These authors recognized that this line, attributed to a forbidden line 
of neutral O. is also blended with Ni I. with the result that the O abundance of 8.83 in Grevesse 


and Sauval ( 19981 ) (8.93 in lAnders and Grevessel . 1 1989( 1 was revised down to 8.69 ±0.05. A number 
of subsequent papers a nalyzing other allowed and forbidden lines in O I, and also molecular OH 


supported this change (lAsplund et all l2004 2005a 


Melendez 


20041: Socas-Navarro and Norton 


20071; Melendez and Asplund], 2008), though Avres 

20081) and Caffau et al. (2008|) offered more 

cautious views. Joining, O, C (Allende Prieto et al. 

2002j; Asplund et all 2005 y; Caffau et al ., 

20101) and N (Oaffau et al. 200911 also underwent downward revisions in their abundances. 

In Table fT| we collect the recommended solar photospheric abundances of various commonly 

observed elements from Grevesse and Sauval (19980. Asplund et al. (2009s). (also given in Grevesse 

et al.. 2010l) Caffau et all (2011), and Scott et al. (20141 

a) and Grevesse et al. (2014|) for comparison 

and reference. Not given here are result from Loclders 

(2010). who for the elements of most interest 

here appears to quote the average of Asolund et al. 

(2009) and Caffau et al. (2011). The data 


sources for each element are all listed in these reviews. For future reference we remark that 
aside from the revisions above, the noble gases will be of interest t o us, since their photosp heric 
abundances a re often determined from coronal observations (e.g., Feldman and WidinsJ . 1990; 
Young, 2005a), on the assumption that no fractionation occurs between the photosphere and the 
corona for these elements. This is a n assu mption we shall scrutinize. Only He is determined 
independently, from helioseismology ( Basu and Antial . 2004 ). in an analysis that includes the 
revised solar metallicity. 


2.2 Helioseismology 

The re-evaluation of the sola r composition above h as revised the s olar m etallicity d own from 0.0170 
iGrevesse and Sa uvall. 19981 to the range 0.0153 ( Caffau et a,L . 201 il l to 0.0134 ( Asplund et al. 


woow cirrei L>uruva|j, ijuiiip a o.uxoo vy aura- u u.o aot: uiiOMiuuu o v u/t/.i , 

2009]), with Asplund et all (2005c) giving a value as low as 0.0122, largely driven by improved 


analysis of lines of C, N, and O. The lower metallicity decreases the sound speed at the base of the 
solar convection zone, yielding now a l arger disagreeme nt between observed and modeled sound 
speeds (see, e.g., Figure 1 in lGuzik and Mussack . 20101 1. The reduced abundances also decrease 
the depth of the convection zone in so la r m ode ls, a gain worsening agreement between models and 
helioseismic inversions (e.g., Basu and Antia, _2004j). 

One early solution proposed was to increase the abundance of Ne to compensate ( Antia and 
Basu, 20051 : Bahcall et all 2005al l. bv bringing the solar metallicity back to its prior value. A 
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Table 1: Recommended Solar Photospheric Abundances of Common Elements 


Element 

Grevesse and Sauval 11998) Asplund et al. 12009) Caffau et al. 12011) 

Scott et al. 12014b.a) 
Grevesse et al. 12014) 

H 

12.00 

12.00 



He 

10.93 ± 0.004 

10.93 ±0.01 



C 

8.52 ±0.06 

8.43 ± 0.05 

8.50 ± 0.06 


N 

7.92 ± 0.06 

7.83 ± 0.05 

7.86 ±0.12 


O 

8.83 ± 0.06 

8.69 ± 0.05 

8.76 ±0.07 


Ne 

8.08 ± 0.06 

7.93 ±0.10 



Na 

6.33 ±0.03 

6.24 ±0.04 


6.21 ±0.04 

Mg 

7.58 ± 0.05 

7.60 ± 0.04 


7.59 ± 0.04 

Al 

6.47 ±0.07 

6.45 ± 0.03 


6.43 ± 0.04 

Si 

7.55 ± 0.05 

7.51 ±0.03 


7.51 ±0.03 

P 

5.45 ± 0.04 

5.41 ±0.03 

5.46 ± 0.04 

5.41 ±0.03 

S 

7.33 ±0.11 

7.12 ±0.03 

7.16 ±0.05 

7.13 ±0.03 

Cl 

5.5 ±0.3 

5.50 ±0.30 



Ar 

6.40 ± 0.06 

6.40 ±0.13 



K 

5.12 ±0.13 

5.03 ±0.09 

5.11 ±0.09 

5.04 ± 0.05 

Ca 

6.36 ±0.02 

6.34 ± 0.04 


6.32 ±0.03 

Ti 

5.02 ± 0.06 

4.95 ±0.05 


4.90 ± 0.04 

Cr 

5.67 ±0.03 

5.64 ± 0.04 


5.62 ±0.04 

Fe 

7.50 ±0.05 

7.50 ±0.04 

7.52 ±0.06 

7.47 ±0.04 

Ni 

6.25 ±0.04 

6.22 ±0.04 


6.20 ±0.04 

Kr 

3.31 ±0.08 

3.25 ±0.06 


3.23 ±0.06 


survey in nearby active star s had previously yielded the abundance ratio Ne/O in the range 0.3- 
0.4 ( Drake and Testal. 2005 b significantly higher than the solar coronal value and closer to that 
suggested. Drake and Testa ( 2005! ) argue that the relative constancy of this ratio in a sample 
of over 20 stars suggests that the “true” Ne/O abundance ratio should be around 0.4, and that 
the variation in the solar corona must be due to some unknown fractionation. More recent work 


sugge sts than an increased Ne abundance (of about 0.5-0.67 dex) is not a complete fix (iLin et al. 
2007|), but a more modest increase of 0.45 dex is still acceptable (see introduction of Guzik and 


Mussack, l2010l) . Ne, and also possibly Ar, are the focus of revisions to composition because having 


in astrophysical sources as proxies for the solar photosphere. 

Asplunc 

et al. 

(20Q3) argue for an 

Ne/O abundance ratio of 0.175 ± 0.031 following the measurements of 

Young 

12005a) in the solar 

corona. Other authors have suggested that increased opacity ( 

Ohristensen-Da 

sgaard et al. 2009!; 


ISerenelli et al . , 2009h , beyond that implied by the increas e due t o data from the Opacity Project 


replacing older OPAL radiative opacities (Bahcall et al., 120051 ) might ease the problem, while 


variations to the solar evolutionary hi s tory with extra episodes of mass loss of accretion h ave also 
been considered ( Guzik and Mussack . 2010l : Serenelli et al . 2011 ). IVillante et al\ ( 20141) provide 
a recent evaluation, accounting for helioseismic and s o lar ne utrino data, again favo uring the older 
photo spheric abundance set of lGrevesse and Sauvall ( 19981 ). Even more recently, IShearer et al\ 
( 20141 ) study the variation of Ne/O measured by Ulysses/SWICS and ACE/SWICS over the solar 
cycle between 1998 and 2012. Their results also favour a low Ne/O abundance ratio in the range 
OJlDOlTS^althoughwithnm expec ted variation. 

Bergemann and Serenelli ( 2014! ) review these and other poten tial solutions to the “solar abun¬ 
dance problem”. Other ideas include that of Lopes and Silk (2013), who revisit the helioseismology 
problem in the light of revised photospheric abundance for the Sun and solar-like stars. Relative 
to solar analogs without planetary systems, the Sun appears to be underabundant in metals, more 


4 















































































so in refractory elements than in volatiles. The mass of the “missing” elements from the solar 
convection zone appears to be similar to the combined mass of the inner terrestrial planets, Mer¬ 
cury, Venus, Earth and Mars, leaving open the possibility that the metallicity of the solar interior, 
specifical l y the r adiati on zone, is higher than that of the convection zone, as appears to be required. 
Lopes and Silk 1 2013 1 explore several such models, computing the solar neutrino spectra with a 


view to future observational capabilities. Their m odels however appe ar to nee d high mass loss from 
the young Sun, which might be problematic (cf. IWoodl . 2004 . 200 0i. Zhanel ( 2014 1 considers the 
effect of a turbulent kinetic flux within the solar convection zone. A negative flux, i.e. turbulence 
propagating from the convection zone to the radiation zone, goes some way to restoring agreement. 
A downwards turbulent kinetic flux requires larger outward radiative and convective energy fluxes, 
which lead to a deeper boundary between the convection and radiation zones. 

Finally the solar helium abundance is most accurately determined from helioseismology, from 
the anomaly produced iimthe sound speed at the depths in the convection zone where helium 


ionizes. 


Villant e et al. (l2014Tl g ive the value as Y = 0.2485 ± 0.0035, apparently the mean of the 


results derived by[Basu and Antia (12004H using GONG or MDI data. 
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3 The Solar FIP Effect: Overview 


As mentioned above, the observation of element abundances in the solar corona and wind has a 


history s 

panning 

Feldman 

, 1992; 


. (IMe; _ 

1992 : Feldman and Laming . 2000l : Feldman and Widinel 20021 . 2003 . 2007h . Here, rather 
than provide an exhaustive review, we attempt to summarize and update the status of solar coronal 
abundances, and refer readers back to these reviews for extensive references. 

In contrast to photospheric abundances which are generally measured and given relative to 
H (i.e absolute abundances), with a few exceptions to be noted below, most coronal abundance 
measurements and results are relative, in that one minor ion is compared to another minor ion 
(often O) with no reference to H. This arises due to the difficulties in observing H. H generally 
has no observable emission lines from the solar corona in bandpasses in which other elements are 
observed, and requires different in situ instruments for detection to those for heavy ions. Both 
difficulties lead to cross calibration issues. Exceptions occur in the UV where H emission lines can 
sometimes be detected, or in the X-ray region where the thermal bremsstrahlung continuum can 
be taken as an indicator of the H abundance. We will highlight these features where appropriate 
below. 

Theoretically, it will also turn out to be easier to discuss relative abundances. Modeling the 
action of the ponderomotive force on minor ions in an ambient H atmosphere is much more straight¬ 
forward than the multifluid calculation that would be necessary to treat the back reaction of the 
Alfven waves on the H fluid itself, and to incorporate those effects into the predicted FIP frac¬ 
tionations. In fact in many cases, it turns out the H and O behave similarly because of their 
strong charge exchange coupling, due to their similar ionization potentials, and this provides some 
justification for our approach. 

In the quiet solar corona and slow speed solar wind, elements with first ionization potential 
(FIP) below about 10 eV are enhanced in abundance by a factor of about 3, with typical variations 
in the range 2-5. This is established by both remotely sensed (i.e., spectroscopic) and in situ 
measurements and at the time of writing, despite early controversy, is now considered an established 
fact. Variations in the fractionation with solar region exist. Coronal holes, and the fast solar wind 
emanating from them are known to have a significantly smaller degr e e of FIP fractionation than 
the quiet corona and slow wind (e.g.. iBochsleri . l2007al : FeldmarJ . 1 998bl ) . iBrooks and Warrenl (|201lh 
show that FIP fractionated abundances measured in a n activ e region with the Extreme Ultraviolet 
Imaging Spectrograph (EIS) on Hinode ([Culhane et an . 2007) match those detected a few days later 
when the s olar wind from this active region would have reached 1 AU. The Sun viewed “as a star” 

( Laming et all Il995l ) sho ws FIP effect only for temperatures above about 10 6 K. This is illustrated 
in Figure [I] taken from Laming et al. (|_1995 ), where the full di s k sola r emission distribution is 
determined from the solar spectrum of iMalinovskv and Heroux (1973), based on emission lines 
from low FIP (filled symbols) and high FIP (open symbols) elements. The two emission measures 
coincide up to logT ~ 6, and thereafter diver ge, indicating overabund a nce o f low FIP e l ements . 
Thi s result is corroborated b y lYounel ( 2005bl ldb iFeldman and Widina (1993); Feldman (1998a), 
and Young and Masonl ( 1998h . and leads to the conclusion that structures emitting at temperature 
below 10 6 K m ust be distinct entities from those responsible for the higher temperature emission. 


Feldma n and Laming! (19941) elaborate this argument based on earlier observations, most notably 


those in Feldmanl ( 19831 . 119871 1. It appears likely that the high temperature FIP fractionated plasma 
resides in coronal loops fed by chromospheric evaporation (t he response of c hr omospheric plasma to 
heat released in the corona conducted downwards, see e.g 


tnc response ot cnromospneric plasma to 
Brav et all . Il991 - L Reale I 2 OI 4 I), whi le the 


lower temperate plasma, attributed to “unresolved fine structu res” bv IFeldman ( 19831. 1987h . can 


now be identified wit h the recently discovered Type II spicules (D e Pontieu et all 1201 ll: Martfnez- 


Sykora et al.. l201ll f. or other dynamic loop structures (jHansteen et all l2014) l recently observed 
with the Interface Region Imaging Spectrograph (IRIS). Individual features distinct from Type II 


6 























































































































spicules at trans ition region temperatures (10 4 K < T <10 6 K) may show FIP fract i onation (see 


exam ples in e.g. Feldman . 1992 : Feldman and LamineJ. 2 000|: Feldman and Widine . 2002, 2003 


2007h . We emphasize again that the result of Laming et al\ ( 1995 1 applies to the “Sun as a star 


Active regions and flares can also have different fractionation, often reduced from that observed 
in the quiet SumThisjsjnost^clearlvseenin^Dectim^cciuire^^fjDlasnMabovej^sunsjao^JlEeldman 


and Widing, 1990t l. [Phillips et al . (Il994lf also measured a maximum Fe enhancement of a factor 
of 2 between photosphere and corona, by comparing the photospheric Fe K /3 line excited by 
flourescence with the collisionally excited coronal Fe XXV res onance line observed during flares 
by the YOHKOH Bragg Crystal Spectrometer ( Culhane et al . 1991 1. Later work with RHESSI 


( Lin et all 12004 1 found much str onger enhancements i n the abundance of Fe relative to H from 


line to continuum measurements ( Phillips et aL . 200fill. Abundances of other elements in flares 


relative to H have been measured by RESIK iSvl wester et al. . 2005 1. FIP fractionations for K 


Ar, Cl, S, Si and Al have been given by Svlwester et al. (l2008h . with more detailed results given in 


subsequent papers ([Svlwester et all 



2013. 12014 1. and given in more detail in 


Table[3]for comparison with model results. IWarren ( 2014 1 finds almost no FIP fractionation in 21 
flares observed w ith the EUV varia bilit y Ex periment (EVE) on the Solar Dynamics Observatory 
(SDO), while Del Zanna and Mason (20141 ) find a typical enhancement in the abundance of Fe 
compared to O and Ne of 3.2 in 9 flares observed by t he Fl at Cr ystal Spectrometer (FCS) o n the 
Solar Maximum Mission (SMM1. iBrooks and Warrenl ( 2012h and Widing and Feldman ( 2008l l find 
similar enhancements. 

Coron al mass ejec ti ons (CMEs) observed in situ often have strong FIP fractionation as shown in 
a survey dReisenf eld et ail l2007t l conducted wit h data the Advanced Composition Exp l orer ( ACE 


Cloeckler et al . TpHlSl l and the Genesis mission ( Burnett et al. . 2003l l. Zurbuchen et al. ( 2004 1 and 


Smith et al . I ( 2001 ) show that frequently these large FIP fractionations are associated with the CME 


flux rope, possibly implying that the flux rope forms in the corona and does not emerge preformed 
from the photosphere (in which case it might be expected to e xhib it ph otospheric abundances). 
CMEs can also exhibit mass dependent fractionation ( Wurz et al. . 2000), most likely indicating 
the role of processes other than FIP fractionation in modifying the elemental composition of the 
solar upper atmosphere. 

We have already mentioned the dif f erence in FIP effect between f as t and slow speed s treams 
observed by Ulysses ( Zurbuchen et al. . 19991 : von Steiger et al. . 2000l l. Lepri et all ( 2013h study 
the evolution of abundances of fast and slow solar wind observed in the ecliptic by the Advanced 
Composition Explorer (ACE). While other solar wind parameters such as charge states, magnetic 
fields and freeze-in temperatures show marked variation, the degree of FIP fractionation (i.e. 
element abundances relative to O) does not. Absolute abundances, measured relative to H, do 
show some variation, with the lowest abundances of He, C, O, Si and Fe being observed in the 
slowest wind at solar minimum. I nteresti ngly, their absolute abundances for O support the revised 
photospheric abundance of O (jAsolund et all 120099 . with the highest mean O/H abundance ratio 
in the slow speed wind at solar maximum being 8.68 (in logarithmic notation), with up per a nd 
lower limits of 8.42 and 8 . 94 res pectively, from their Table 1. Previous studies of this sort dBochsler . 
2007bl : von Steiger et al . 2010 1 have favoured the “older” O abundance (8.83 lGrevesse and Sauval 


von Steiger et al. (1201 0) assume that the high speed solar wind gives the most faithful 


representation o f pho to spheric abundance s (their O/H value in slow speed solar wind agrees with 


lAsplund et al . 2009 ;. Bochsler ( 2007b ) plots O/H against He/H in his Figure 3. and shows 


a striking correlation between them, arguing that inefficient Coulomb drag is the cause of both 
variations. This might appear to be proble matic, in that O is not always depleted relative to H 
(unless one takes oldest O/H value from 1 Anders and G re vessel . 1 198()l . as the correct value), but He 
is, with the problem becoming worse for the more recent O/H values. 

The solar slow speed wind helium abundance has long been known to be depleted from the 
photospheric value, and this depletion is now established to vary with solar slow wind speed and 
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Figure 1: Full disk solar emission measure distribution, with low FIP elements depicted by solid 
symbols, high FIP elements by open symbols. The transition from photospheric abundances at 
logT < 6.0 to coronal abundances at higher temperatures, indicating that di fferent types o f coronal 
structures are emitting above and below this temperature. 

Reproduced by permission of the AAS. 


Figure from Laming et al. ([19951 1. 
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Figure 2: Left: Helium abundance relative to hydrogen in the slow speed solar wind over 1.5 solar 
cycles. The curve colors denoted the wind speed at which at abundance is measured. He is more 
variable in the slowest slo w speed solar w ind. The black curve gives the monthly smoothed sunspot 
number. Figure from Kasper et alj (J2Q12). Reproduced by permission of the AAS. Right: Helium 
abundance relative to oxygen in the slow speed solar wind, measured by ACE/SWICS. The same 
trend of depletion with wi nd speed as for He/H is se en, but the solar cycle dependence is less 
pronounced. Figure from Rakowski and Lamind ( 2012 1. Reproduced by permission of the AAS. 



Intensity 



Figure 3: Maps of intensity, non-thermal velocity, photospheric magnetic field and coronal FIP bias 
derived for a sigmoidal anemone active region observed by EIS/Hinode. The overlaid black ellipses 
show the footpoints of loops, coincident with strong photospheric magnetic field concentrations, 
while the dotted black lines show the loop connections between them. The strongest FIP effect 
can b e seen at the loop footpoints, also coincident with strong nonthermal velocities. Figure from 
Baker et al. ( 2013l l. Reproduced by permission of the AAS. 
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the phase of the solar activity cycl e ( Aellig et all 1200 ll : iKasner et all 2007 ) in a s i milar manner to 
that mentioned above ( Lepri et aill2013l ). Figure 121 (left panel. from lKasper et al. . 2012 1 shows the 
He/H abundance ratio measured by the Wind spacecraft over 1.5 solar cycles, in bins of different 
solar wind speed. Greater variability is seen in the slowest wind speed bins, and in the slowest 
bin, the variability matches the smoo the d sunspot number. The right panel of Figure [2] show a 
similar study with ACE bv iRakowski and Laming! ( 2012 1. where the abundance ratio He/O is 
plotted. Broadly similar behavior with solar wind speed is seen, but the solar cycle dependence 
is less pronounced, but still apparent. The fast solar wind helium abundance is also depleted, 
but to a lesser degree than t he s low win d, and is also less variable. Measurements of the helium 
abunda nce in co ronal holes (Laming and Feldmanl. 2003ll and in the quiet solar coron a (Laming and 
Feldman, 200il l generally show similar values to those measured in the solar wind, indicating that 
the helium depletion occurs l ower do wn in t he so lar atmosphere than the corona. A spectroscopic 
measurement in a solar flare ( Feldman et al. . 2005l l showed a much higher abundance of helium, in 
qualitative agreement with the by now routine observation in situ of en hanced (i.e., less depleted) 
helium abundance in CMEs (e .g., IWimmer-Schweingruber et al\ . l2006ll . 

Widing and Feldman ( 200ll l studied the variation of the FIP effect with time in newly emerged 
active regions. They found that the new loops emerged with photospheric abundances, and grad¬ 
ually changed to coronal abundances, i.e., developed a FIP effect, over the course of a few days. 
Consequently a new active region should hav e a weak F IP effec t, while an older one would show 
strong fractionation. Recent observations by Bakey et al} (120131 1 corroborate this view. Figure [3] 
shows maps of intensity, non-thermal velocity, photospheric magnetic field and coronal FIP bias 
derived for a sigmoidal anemone active region observed by EIS/Hinode. The overlaid black ellipses 
show the footpoints of loops, coincident with strong photospheric magnetic field concentrations, 
while the dotted black lines show the loop connections between them. The strongest FIP effect can 
be seen at the loop footpoints, also coincident with strong nonthermal velocities. Weaker FIP effect 
is seen along the loop connections, suggesting that the FIP effect originates in the chromosphere, 
and is communicated to the coronal loop by transport processes. The active region seen here would 
then be considered relatively new, since only weak FIP effect is seen in its coronal connections. 
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4 Stellar FIP and Inverse FIP Effects 


The measurement of element abundances in stellar coronae became possible for the first time 
with the 1992 June 7 launch of the Extreme Ultraviolet Explorer (EUVE) satellite. The use 
of grazing incidence gratings allowed strong lines in the EUV spectra of stars to be resolved, 
allowing the ac quisitio n of dat a on stellar corona of similar quality to the early solar spectra 
analyzed by e.g. JPottaschl (1 9631 ). The Advanced Satellite for Cosmology and Astrophysics (ASCA) 
launched on 1993 February 20 provided X-ray spectroscopy of stellar coronae with CCD-level 
spectral resolution. Individual lines could not be resolved, but the He- and H-like line complexes 
of different elements could. Other missions from which a few results emerged were Ginga (launched 
on 1987 February 5), ROSAT (launched 1990 June 1) and BeppoSAX (launched 1996 April 30) 
which carried proportional counters for spectroscopy (gas scintillation proportional counters in the 
cases of Ginga and BeppoSAX) which allowed measurements of the He-like Fe complex of lines 
with respect to the surrounding continuum. More recently Suzaku (launch 2005 July 10) has also 
provided CCD X- ray spectra, although the higher resolution calorimeter failed soon after launch. 

The review of lFeldman and Laming ( 2000l l covered the status of this new field at a time when 
the results from EUVE and ASCA were available, but before the launches of Chandra and XMM- 
Newton (launched 1999 July 23 and 1999 December 10, respectively). These satellites both carry 
grating instruments as well as CCD cameras, and so high spectral resolution X-ray spectra of stellar 
corona e are acquired with significantly higher throughput than was the case with EUVE or ASCA. 
Testa ( 2010h provided a more updated view of t he field, encom pass ing a wider variety of stellar 


ta rgets. Here we will r ecap t he discussion of lFeld man an d Laming ( 2000h . and the later reviews 
of Favata and Micelal (2003), Drakel (2003) and Testal ( 2010 1. concentrating on stars of varying 
spectral type but in other respects similar to the Sun, with a view to the theoretical discussions to 
follow in Sections [S] and [7] 
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Table 2: List of FIP Bias Assessments for Main Sequence Stars (updated from (Wood et al. 20125). 


Star 

Alternate 

Name 

Spectral 

Type 

Radius 

(Re) 

log Lx 
(ergs s 1 ) 

log PA 

ergs cm -2 s -1 

Ne/Fe 

-^bias 

Ref. 

7r 3 Ori 


F6V 

2.05 

28.96 

5.55 

2.00 

-0.55 

Wood and Laming (2013) 

/3 Com 

HD 114710 

GO V 

1.08 

28.21 

5.36 

0.73 

-0.668 

Telleschi et al. (2005i 

tt 1 UMa 

HD 72905 

G1 V 

0.91 

28.97 

6.27 

1.60 

-0.645 

Telleschi et al. (2005) 

X 1 Ori 

HD 39587 

G1 V 

0.98 

28.99 

6.22 

1.91 

-0.555 

Telleschi et al. (2005) 

Sun 


G2 V 

1.00 

27.35 

4.57 

1.02 

-0.600 

Feldman and Laming (2000) 

a Cen A 

HD 128620 

G2 V 

1.22 

27.00 

4.04 

1.41 

-0.410 

Raassen et al. (2003b) 

k Ceti 

HD 20630 

G5 V 

0.98 

28.79 

6.02 

2.51 

-0.462 

Telleschi et al. (2005) 

£ Boo A 

HD 131156A 

G8 V 

0.83 

28.86 

6.24 

3.80 

-0.344 

Wood and Linskv (2010) 

70 Oph A 

HD 165341A 

K0 V 

0.85 

28.27 

5.63 

2.40 

-0.403 

Wood and Linskv (2006) 

36 Oph A 

HD 155886 

K1 V 

0.69 

28.10 

5.64 

4.90 

-0.250 

Wood and Linskv (2006) 

36 Oph B 

HD 155885 

K1 V 

0.59 

27.96 

5.63 

3.80 

-0.328 

Wood and Linskv (2006) 

q Cen B 

HD 128621 

K1 V 

0.86 

27.60 

4.95 

1.44 

-0.478 

Raassen et al. (2003b) 

e Eri 

HD 22049 

K2 V 

0.78 

28.32 

5.75 

4.68 

-0.050 

Wood and Linskv (2006) 

£ Boo B 

HD 131156B 

K4 V 

0.61 

27.97 

5.62 

7.41 

-0.185 

Wood and Linskv (2010) 

70 Oph B 

HD 165341B 

K5 V 

0.66 

28.09 

5.67 

8.71 

0.138 

Wood and Linskv (2006) 

GJ 338AB 


M0 V+M0 V 

0.59+0.58 

27.92 

5.30 


0.305 

Wood et al. (2012) 

EQ Peg A 

GJ 896A 

M3.5 V 

0.35 

28.71 

6.84 

16.12 

0.450 

Liefke et al. (2008) 

EV Lac 

GJ 873 

M3.5 V 

0.30 

28.99 

7.25 

13.92 

0.474 

Liefke et al. (2008) 

EQ Peg B 

GJ 896B 

M4.5 V 

0.25 

27.89 

6.31 

14.67 

0.417 

Liefke et al. (2008) 

AD Leo 

GJ 388 

M4.5 V 

0.38 

28.80 

6.86 

18.25 

0.536 

Liefke et al. (2008) 

Proxima Cen 

GJ 551 

M5.5 V 

0.15 

27.22 

6.08 

11.98 

0.471 

Liefke et al. (2008) 

Very Active Star Sample (log L 

x > 29) 







EK Dra 

HD 129333 

G1.5 V 

0.89 

29.93 

7.25 

3.71 

-0.277 

Telleschi et al. (2005) 

AB Dor 

HD 36705 

K0 V 

0.79 

30.06 

7.48 

17.11 

0.488 

Giidel et al. (2001a) 

AU Mic 

HD 197481 

Ml V 

0.61 

29.62 

7.27 

24.30 

0.695 

Liefke et al. (2008) 


































Wood and Linsky (2010) conducted a survey of stellar FIP effects restricted to stars with X-ray 
luminosities less than 10 29 erg s _1 , with results shown in Figure Q] (adapted from Figure 9 of their 
paper, also including points for aCen B, 7r 3 Ori and GJ338). By excluding the most active stars, 
a clear trend of decreasing stellar FIP effect with later spectral type is uncovered. The Sun is 
at the bottom left of Figure [2 with a logarithmic FIP bias —0.6 (expressed as log(X/H) phot — 
log {X/H) COI , so that this is a coronal low FIP enhancement of IQ 0 6 ~ 4) , along with other stars 
of similar spectral type (7r 1 UMa (G1V) and y 1 Ori (GOV); iGiidel et oil 2002). The degree of 
fractionation diminishes as one moves to later spectral type, becoming zero at about K5, and an 
inverse FIP effect is observed in the M stars. The magnetic fields of a sunspot umbra and penumbra 
are likely to be similar to those found m ore ubiquitously in the atmospheres of later type (i.e. M) 
stars (seelDonati and Landstreetl. 20091). so the results of Feldman and Widingl (1990h and Phillips 


et al. (119941 ) appear consistent with these stellar results. At the extreme left hand side, the result 
for 7T 3 Ori suggests that the FIP effect saturates at a value similar to that found in the Sun, a low 
FIP enhancement of about 4, and does not continue increasing at spect ral ty pes ea rlier than this. 
The properties of the various stars are summarize d in Table 0 IWood et q.Z. 1 (I2012h show that the 
sample of T Tauri stars from IGiidel et al. ( 2007 1 show a similar trend of Ne/Fe as in Figure U 


but generally at higher FIP bias (i.e., stronger Inverse FIP effect). Other more active stars (X-ray 
luminosities larger than 10 29 erg s _1 ) also show generally larger inverse FIP than those in Figured 
We do not consider these further because of complications to the stellar physics introduced by fast 
rotation and tidal interactions in close binaries that are typical of these types of stars. 

Anticipating a connection between the coronal abundance anomalies, and the chromospheric 
wave Held, arising either from waves produced in the corona and propagating down, or as a result 
of hclioseismic or asteroseismic p-modes, we make the following comments about some of the 
important stars in Figure [4] 

_ a Cen AB: This inactive G2 V+Kl V binary is not in the original sample of lWood and Linskv 

12010 l l. but it is of prime interest since its fundamental parameters (photospheric abundances, p- 


mode frequencies, mass, radius, surface gravity, and effective temperatures) are very well known 
for both stars (e.g.. Bedding et al.. 2004; Butler et a(l 2004; Kieldsen et all 20051: Porto de Mello 

et al.. 20081: Bruntt et all. 2010): 

Karoff et al. 20071: Chaplin et al.. 2009). Coronal abundances were 

first measured by 
two in the unreso 
coronal source, wl 
and Chandra (e.g 

Drake et al. 

ved binary. S 
iich places it 
. Raassen et 

19971) from EUVE data, yielding a FIP effect of about a factor of 
ubsequent observations determined that a Cen B is the dominant 
n its spot in Figure 1. There are many observations by both XMM 
ali 2003b; Licfke and Schmitt. 2006^ and model chromospheres of 


both stars have been developed (|Avres et all 1 197f: : Jordan et o/.Tl987 : Vievtes et al\ . POOfil ). Fe 


XII 1242 A linewidths are ava ilable in lAvres et T7/T )(l.'ji ;. 


e Eri: lLaming et all 1 19961 ) first measured coronal abundances for this moderately active K2 V 
star with EUVE. More recent spectra from XMM and Chandra have been exte n sively analyzed as 
well (ISanz-Forcada et al.. . 20041 : Wood and Linskv . 20061 : Ness and Jordanl . 2008). iDrake and Smitbl 
( lOO.lh give fundamental parameters (see also T a.ble 2 inlVievtes et all 2009). Only theoretical 
estimates of e Eri’s p-mode spectra are available (iGai et al. 1200811 . Several so u rces provide mode l 
chromospheres ( Jordan et al. . 1987 : ISim and Jordanl . 2005^ Vievtes et al . 20091 ). lAvres etal. ( 20031 ) 
give the Fe XII 1242 A width from which coronal turbulence may be deduced. 

£ Boo A: Thisjs_aG8_Vdwarf with coronal abundances first measured with EUVE by Laming 


and D rake (j 19991 ) and IDrake and Kashvanl " 200ll). and later wit h Ch andra (Wood and Linskv , 
2010t ). Model chromospheres are available from lKelch et al. ( 1979t ) and Llordan et ali dl987h . 

70 Oph A : This K0 V dwarf with an X-ray spectrum was anal yzed by IWood and Linskvl 
( 2010 l ). Its fund amental parameters are given by Bruntt et al\ (I201C), inc luding p-mode oscilla¬ 


tion frequencies (Carrier and Eggenberger, 2006; Eggenberger et all 


20081 ). No empirical model 


chromosj]herejs i ^w|ilable i fo]Mhis ii star J 2Du^herejsjrUieoreticalmre_from i Schmitz a nd Ulmschnei- 
der (il980h . Further guidance may come from the K dwarf model chromospheres of Vievtes et al 
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Figure 4: Survey of FIP fr acti o nation (log (X/iL) phot — log ( X/H ) COI ) observed in a sample of 
dwarf stars, updated from Wood and Linskvl ( 2010 1. where values below zero indicate a solar-like 
FIP effect and v alues above zero in dicat e an inverse FIP eff e ct. P oints for GJ338 and n 3 Ori have 
been added f rom Wood et al. ( 2012 1 and Wood and Laming! 1 201 .'ill . respecti v ely, a nd aCen B from 


Drake^t_a(J (jj_997|). Diamonds indicate m easurements from Telleschi_el_aL ( 2005l l. triangles from 


Liefke et al\ ([2008) , with a solar value from Feldman and Lamina! ( 2000f l. For all GK stars the FIP 


bias c alculations include corrections for stellar photospheric abundances from lAllende Prieto et al 
(|20 04), but for the M stars there are in general no stellar photospheric measurements available 
so we have to simply assume solar photospheric abundance apply. Among these, EV Lac is the 
strongest cas e, due to chro mospheric ev aporation of photospheric abundance material observed 
during flares ( Laming and Hwanel . 20091 . and references therein). For the purposes of this figure, 
we avoid extremes of stellar activity, confining our attention to stars with log Lx < 29. Reproduced 
by permission of the A AS. 


( 2009 1. It is also included in the survey of forbidden line observations bv lAvres et al\ (l2003h 
M Dwarfs: The cluster of M dwarfs at top right in Figure 1 are studied by iLiefke et a,l 


{ 2008h . A model chromosphere for AD Leo is gi ven by Fuhrmeister et aj\ ( 2005h , and a grid of 
model chromospheres for Ml dwarfs is given bv iHoudebine and Stempela ( 1997 b lAvres et al 


( 2003l l give linewidths for the Fe XII 1242 A and Fe XXI 1354 A forbidden lines. EV Lac is 
most recently studied by Laming and Hwang] ( 2009f ) using Suzaku observations before and during 
a flare. It exhibits abundance change during flares; the inverse FIP quiescent corona giving way 
to a more “normal” (i.e., solar photospheric) abundance pattern. This is interpreted as evidence 
of chromospheric evaporation, where unfractionated plasma is evaporated up into a flaring loop, 
lending more confidenc e to th e inv erse FIP interpretation. 

Favata and Micela ( 2003h and Sanz-Forcada et~ai\ ( 2004h caution that for several stars with 
apparent metal depletion in their coronae, determined from the Fe/H ratio, the coronal abund ances 


mere ly reflect metal poor photospheres. In Figure IH AD Leo maybe such a case ( .Tones et al. 

11996 1. although it clearly shows a nonsolar Ne/Fe abundance ratio. However metal depleted 
or inverse FIP effect, do appear to clearly exist in some cases, e.g., II Peg (K2 IV plus 


coronae, 
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an un seen companion: iHuenemoerder et al. 120011 ). AR Lac (G and K subgiants in a 1.98 day 
orbit; Huenemoerder et all 1 200. 'll) and AB Dor CK2 IV-V with a 0.515 day spin period: Sanz- 


Forcada et al., 2003i ). Further, the variation of element abundance during stellar flares, in "which 
initially metal depleted plasma evolves towards the standard composition, interpreted in terms of 
the chromospheric evaporation of unfractionated plasma to the coronal flare site, seems to require 
the exi stence of such abund ance anomalies. Such phe nomena are observed in HR 1099 (K1 IV and 
G5 IV; Audard et, a/1 2001 1 Algo l (B8 V and K2 IV; Favata and Schmitt . 1999t h and UX Ari (G5 
V and K0 IV; Giidel et al. I . 19991 ) w here the later type s ubgiant is taken to be the main source of 
coron al emission, as well a s AB Do r (iGiidel et all 2001b! ) . YY Gem (dMe and dMe with 0.814 day - 
orbit; Giidel et aLl2001al h II P eg ( Mewe et al.. 1997), AT Mic (dM4.5 and dM4.5; Raassen et al. 


2003al ). and EV Lac (dMe3.5e: ILaming and Hwangi . 12009 1. Many of these more active stars show 


stronger Inverse FIP effects than shown in Figured] 

Not included in Figure dl but also of interest for having neither a solar-like FIP effect nor an 
inverse FIP effect is P rocyon. Coronal abundances measured with E UVE, Chandra and XMM 


( Drake et al. . 1995 b a: Ra assen et al. I . 2002 : Sanz-Forcada et al . 2004 1. At spectral type F5 does 


not follow the trend of Figure 1 (it is, however, a subgiant, not a solardike dw arf st ar ). Extensive 
asteroseismological observations are available for Procyon (e.g.. iMosser et al. . 20081 : Leccia et al. . 
[2007 ). In fact, its p- mode lifetimes are known to be significantly shorter than those of the Sun 


(Bedding et al.. 2010h . Chromospheric models are available for u s to use for this well-studied star 
( Avres et al . 1974 : Evans et al . 1975 : Brown and Jordanl. 198 il l. 


Another interesting case is r Bootis (Maggio_et_aL 120111 1. It is a F7 V dwarf, with a M2 V 


companion. Following the assumption of Maggio et all ( 201ll ) that the F7 V star is responsible for 


the X-ray emission, it is also discrepant from the trend in Figured! with Fbias = —0.17 at spectral 
type F7. r Bootis A has already attract ed interest be cause it hosts a close-in giant planet (r Boo 
b: Prot = 3.31 days, Msini = 3.9Mr) feutler et~al. 


19971 ). A star-planet interaction has been 


noted, in that the plan et appears to b e inducing an active region on the star, that leads the planet 
by ~ 70° in longitude (Sh kolnik et all [2008 : Walker et al. . 2008 ). and this opens the possibility of 
the planet also affecting the coronal abundances, due to waves induced in the stellar atmosphere. 
Another possibility of course is that the XMM observation is contaminated by emission from the 
M2 companion. Another c uriosity is that t he ph otospheric abundances of r Boo A do not agree 
with the trend repor ted bv lLopes and Si lk (2013), being of order 0.2 dex higher than solar values 


( Maggio et al\ . 2011). e Eridani also hosts a Jupiter mass planet, but shows no coronal abundance 


discrepancy from the trend in Figured! 
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5 Early Theoretical Models: Overview 

5.1 Diffusion Models and Variations 


The earliest attempts to explain that FIP effect invoked various processes like thermal and ambipo- 
lar diffusion (in a stationary atmosphere) or inefficient Coulomb drag (in a moving case). Thermal 
diffusion arises in a temperature gradient, with minor ions diffusing towards higher temperatures. 
This happens because collision cross sections decrease with increasing collision energy, and so the 
force from the hotter particles from the direction in which the temperature is rising is lower than 
than from the colder direction. This obviously will select ions over neut rals, and accelerate them 
towards the corona. However such a process is inherently slow (IHenoux . 1995 L 1998), and requires 
static conditions for periods ra nging from tens of hours ( Hansteen et all 1994 1 to days or weeks 
(Killie and Lie-Svendsen, 120071 ). This appears increasingly at odds with the modern view of the so- 
lar chromosphere as a dyna mic environment, continually being pertur bed by the passage of shocks 


(|Carlsson and Steinl . 12002) or reconnection at chromospheric layers (llsobe et all 12008). It also 


conflicts with the argument that in the absence of mass supply from below, the solar wind would 
empty the corona within 1-2 days. Clearly, any fractionation mechanism that takes longer than 
this to change coronal abundances cannot be right. Ambipolar diffusion refers to the diffusion of 
neutrals along an ionization balance gradient, and must also b e con si dered i n such models, though 
by itself does not appear capable of causing FIP fractionation. lHansteen et al. (1997) also demon¬ 
strate that chromospheric mixing, in their case between hydrogen and helium, is also necessary to 
prevent gravitational settling of helium in the chromosphere and to give a realistic abundance of 
helium in the solar wind. 

In an attempt to speed up the fractionation process by thermal means, various authors have 
modeled the separation of elements entrained in a flow of neutral hydrogen and protons, von Steiger 
and Geiss (| 19891 ) considered such a flow driven across magnetic field lines, either a horizontal flow 
across vertical field lines, where gravity plays no role in the fractionation, or a vertical flow driven 
across horizontal field lines, where it does. FIP fractionations matching observations reasonably 
well are achieved, but the models themselves must be considered highly idealized and unlikely to 
represent the real sun. No mechanism is suggested to produce such a cross-fiel d flow . In a tre atment 
of fractionation in a (more plausible) vertical flow along vertical field lines, iMarsch et all ( 1995 ) 
solve diffusion equations for ions and neutrals separately in the background flow. These equations 
only include photoionization of neutrals. No recombination of ions and electrons is accounted 
for (their equations 6 and 9). At the lower boundary the gas is assumed completely neutral, 
with density n (s = 0) = no- The neutral density gradient is specified at the upper boundary, 
taken to be where the gas is completely ionized, and is accelerated into the solar wind. Here 
dn/ds (s = S) = 0 and the ion density also n + (s = S) = 0. Finally the ion density gradient at the 
lower boundary dn + /ds (s = 0) = 0. The boundary conditions on n and n + are justifiable, and do 
not greatly affect the solutions. Different choices can lead to the same diffusive fluxes. The choices 
of boundary conditions for the density gradients are less clear, and these are crucial to obtaining 
FIP fractionation in the model. In fact it is precisely the choice of dn + /ds (s = 0) = 0 that 
specifies the ion flux at the upper boundary, with value no\/ D/t where D = v?/ vjh, the diffusion 
coefficient for neutrals of element j in background of neutral H, in terms of its thermal speed 
and collision frequency with H, and r is photoionization time for neutral j in the chromospheric 
radiation field. The fact that the ion flux at the top of the chromosphere depends ultimately on 
the characteristics of neutrals at the bottom through the diffusion equation, and especially so for 
the FIP fractionated low FIP elements that in reali ty have very sma ll neutr al fra ctions throughout 
the region of interest, is indicative of problems. McKenzie et al\ ( 19981 ) and McKenzi e ( 2000t ) 
point out that in a gas where collisions play the dominant role in coupling species together, FIP 
fractionation cannot occur in a one dimensional steady state model. If all elements enter at the 
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lower boundary with the speed of H, and leave at the upper boundary with at the proton speed, 
continuity demands that no fractionation occur. They also conclude that the boundary condition 
on the density gradient at the lower boundary is the cause of the fractionation in this model, and 
that su ch fra ctiona tion must have taken place below the lower boundary. Indeed, equation 71 of 
Marsch et al. ( 1995h giving the upper boundary condition on the exiting ion fluxes suggests that the 


FIP fractionation is already embodied in the inputs to the model, rendering the model explanation 
spurious. 

Fractionation in such a class of models is not completely ruled out. If the chromosphere can 
respond to the enhanced coupling of ions to protons in the upper chromosphere (compared to that 
of neutrals) by allowing diffusion to supply ions from below at a greater rate, i.e. dn/ds has to 
vary with species. But this reliance on diffusion places a limit on the spee d with which abundance 
modification may take place. There are some newer models ( Pucci et al. . 201 0) investigating the 
changes in abundances with hydrogen flux through _a chromospheric magnetic funnel, but this 
requires significant fine tuning for each element ( Pucci et al. . 2010l . consider O, Ne, and Fe), and 
appears unlikely to be able to get every element right at the same time, even if the a tmosphere 
was sufficiently quiescent (see also Bvhring et al. . 2011 : Bvhrinel . 201 il l. B &et al\ ( 201, il l make the 
point that above the chromospheric temperature minimum, the atmosphere is convectively stable, 
but ignore the fact that waves from the convect ively unsta ble re gions may continue to propagate 
upwards to perturb higher altitudes, as in e.g. Heggland et al\ ( 201ll ). They consider the effect 
of gravitational settling within the chromosphere as a means of fractionation among O, Ne, S, 
and Fe. The mechanism only appears capable of producing a depletion of high FIP elements (no 
enhancement of low FIPs as observed), and then only with a static atmosphere. An episodic upflow 
to supp ly t he cor ona, pu nctu ated by periods of stasis to allow the fractionation to occur. Both 
IPucci et ail ( 2010 h and B det al\ ( 20131 ) use somewhat unrealistic chrom ospheric models, in that th e 
degree of hydrogen ionization for a given density is lower than that in Avrett and Loeseil ( 20081 ). 
at l east f or the models which come closest to matching observed abundance anomalies. Although 
the I Avrett and Loeser (2008) model is quasi-static designed to match s pectroscopic observati ons, 
it matches quite well with average chromos pheric structures modele d inlHeggland et al . I d201lj) , as 
does its antecedent ( Vernazza et al . 198ll . VALC) compared with Carlsson and Stein (l2002h . at 
least in terms of electron density. 


The helium depletion explained by such processes ( Bvhrin J . 2011 ) is also at odds with the 
observ ation that it (and the o t her minor ions) flow faster than hydrogen in the solar wind at 
1 AU ( Neugebauer et a,l. . 199fil : Kasper et al. . 20081 : Bourouaine et all 2011 ). A further s tage of 
solar wind acceleration must set i n at a level a bove that where the He depletion sets in ( Wand . 
20081 ). Noci et al. 1 19971 ). following I Geiss et~al. 1 1970l ). make the point that if such processes were 


the origin of the solar wind He abunda nce depletion, other heavy ions should also be depleted, 
which is now known not to be the case. iBochsler ( 2007b ) revisits this, comparing He, O, and Ne 
for which He appears most affected by gravitational settling, or inefficient Coulomb drag, in the 
cor ona . In th is case the He abundance should correlate with the proton flux, which is not observed 
((Wang], 120081 ). However all these elements appear also to vary with respect to H, leaving open 
the possibility that quasi-thermal or diffusive processes play some role in establishing the absolute 
coronal abundances. 


5.2 Thermoelectric Driving 

Perhaps the best early model was that of Antiochosl ( 1994 ). Cross B diffusion of chromospheric 
ions into a flux tube by a thermoelectric force associated with downward heat conducting electrons 
from a closed loop enhances the loop footpoint abundances of ions, but not neutrals. Given 
V x E = 0 in steady state conditions, the transverse gradient of the longitudinal electric field (due 
to the plasma resistivity) must be balanced by a longitudinal gradient of the transverse electric 
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field. This transverse electric field points into the flux tube, thus concentrating ions therein. The 
predicted fractionation pattern is proportional to ion mass 1 / 2 (it is independent of mass for the 
ponderonrotive force). In a sense, this is conceptually very similar to the ponderonrotive force 
model, in that it is the loop’s response to coronal heating that causes the fractionation, except 
here it is mediated by heat conduction rather than Alfven waves. The main problem is that heat 
must be conducted down to regions of the chromosphere without transverse magnetic structuring 
(there must be plasma surrounding the flux tube), and this appears unlikely. It is unclear where 
in the chromosphere horizontal structuring sets in, but a natural place to expect it is likely to 
be the equipartition layer, where the Alfven speed and sound speed are equal, roughly where 
the plasma /3 ~ 1. The stopping distance of 1 keV electrons is of order 1 km at a density of 
10 12 cm -3 , so downwards heat conduction is unlikely t o reach this far, and o t her authors have 
recen tly considered energy transport by Alfven waves ( Fletcher and Hudson! . 20081 : iHaerendel 
2009). Alfven waves may still fractionate material in the upper magnetically structured layers of 
th e chromosphere, b ecause ions move vertically along the flux tube rat her than horizontally across 
it. Antiochoa (1994) also does not allow for an “Inverse FIP Effect”. iLaming and Hwangl ( 20091 ) 
discuss how the ponderomotive force may also inhibit downwards heat conduction, but large Alfven 
wave amplitudes are required. 


5.3 Chromospheric Reconnection 

Arge and Mullan dl998f l consider chromospheric reconnection giving ions a larger density scale 
height than neutrals. They model this using the Zeus code, iterating between its hydrodynamic 
and magnetohydrodynamic implementations to treat the neutrals and ions respectively, with these 
two fluids being coupled by collisions. They only calculate Si and Ne as examples, so it is not 
possible to evaluate the full fractionation pattern such a process would produce. The ion-neutral 
coupling rate they use appears to be that due to charge exchange collisions between protons and 
neutral hydrogen (their equation 6), and this appears to be applied to all elements, with possible 
modifications to the collision velocities. This is_not correct, an d the coupling should more properly 
be described by elastic scattering (cf. iMalvshkin and Zweibeh 20111 ). In fact, looking at these ion- 
neutral collision rates a mass depe ndent fractionation appears likely, n ot the FIP effect , whi ch might 
be re levant to the observations of IWurz eiToZ ( 2000f ) 3 By contrast lLamin 5 (l2004al f and [l arningi 
(2012) devoted considerable effort to the correct description of these atomic processes. Arge and 


Mullan ( 19981 ) also predict FIP effect everywhere, with no difference between coronal holes and 
closed loops, unless extra assumptions are made about the presence or absence of chromospheric 
reconnection in open and closed field, and cannot reproduce an “Inverse FIP Effect”. 


5.4 Ion Cyclotron Wave Heating 

Schwadron et al. ( 1999 1 offered the first mechanism of fractionation by wave-particle interactions. 
Coronal ion cyclotron (IC) waves propagate down and heat chromospheric ions, not neutrals. The 
FIP effect arises from the preferential heating combined with the coupling of the chromospheric 
minor ions to a background flow of H atoms and protons. The formalism for this coupling is elegant, 
and forms the basis for models invoking the ponderomotive force to be discussed at greater length 
elsewhere in this paper. 

The fractionation pattern depends on the assumed spectrum of IC waves, since a resonant 
wave-particle interaction is assumed. This is not specified in their paper, except through the wave 
interaction rate given as u sw = N w Cl s where f l s is the gyrofrequency of ions of element s, and N w 
is a constant of order 70. Equating this to the pitch diffusion coefficient, D = (n/A) fi s (SB 2 /1? 2 ), 

1 This is similar to the case in lDrake et all : 200!lh . where pick-up ion behavior leads to mass dependent fraction¬ 

ation. 
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where 5B is the magnetic field perturbation due to the wave, and B is the ambient magnetic 
field, we find SB/B ~ loH With B ~ 10G and a density of order 10 10 cm 3 , the wave velocity 
amplitude is 2000 km s _1 . Of course these estimates should be treated with caution, because they 
are derived by applying quasi-linear theory well beyond its regime of validity, but the essential 
point that the wave energy requirements of this model are implausible remains. If the coronal ion 
cyclotron wave derived form a turbulent cascade, then even higher wave amplitudes should be found 
at lower frequencies, assuming for example a Kolmogorov spectrum. Loop resonant frequencies 
are lower than ion cyclotron frequencies by several orders of magnitude (typically 6-8), leading to 
a huge increase in the intensity of these waves if a -5/3 spectrum is assumed. Further, waves in 
the ioncjdotronGreauenc^^ange^ayidl^dam^^^charge^xchange^eaction^ejg^Kulsrud and 
Pearce, 19691 ). with rate 1 - 100 s -1 , depending on the neutral fraction. Thus waves traveling 
at the chromospheric Alfven speed are damped after traveling a few km, increasing the energy 
requirements still further. 


5.5 Summary 


To summarize, while some of the models described above might come close to describing the solar 
FIP effect, they do so for rather contrived magnetic field ge omet ries or for various other extreme 
assumptions. In all cases except those of Antiochosl 1 1994 ) and Schwadron et all (1999), whose 
mechanisms both derive from the coronal response to energy deposition, no natural account of the 
difference between closed and open field is given. Further, none of the models appear capable of 
explaining the Inverse FIP Effect._ 

Given the relative success of Antiochosl ( 1994h and I Schwadron etal. ( 19991 1. it seems reasonable 
to pursue the byproducts of coronal heating as the key to the fractionation. Given the high 
frequency waves do not work due their energy requirements, and that lower frequency Alfven waves 
carry far more of the wave energy, the next question to ask is how such waves could interact with 
chromospheric ions. Clearly a resonant interaction is not possible, but such waves can interact 
nonresonantly through the ponderomotive force. This is discussed more fully in the following 
sections. 


2 Bohm diffusion corresponding to N w ~ 1 would give 5B/B ~ 1. 
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6 


The Ponderomotive Force Model 


6.1 Overview 


We now turn to a description of our model of the FIP effect, which invokes the ponderomotive 
force due to magneto-hydrodynamic waves as the agent that separates io ns from neutr als. Pon¬ 
deromotive forces in magnetospheric and space plasmas are reviewed by iLundin and Guglielmi 
(2006). According to them, “Ponderomotive forces are time-averaged nonlinear forces acting on 
media in the presence of oscillating electromagnetic fields. The word ponderomotive comes from 
the Latin words pondus (ponderis) meaning “ heaviness ” and motor.” The complex dynamics of a 
system acting under Lorentz forces may be considerably simplified by averaging over the period 
of oscillations and described instead by ponderomotive forces. They can be seen to arise from the 
effects of wave refraction in an inhomogeneous plasma. In a nonmagnetic plasma, the refractive 
index, is given by e = 1 — w 2 /u; 2 where co p is the plasma frequency. Waves are refracted to high 
refractive index, which means low plasma density. The increased wave pressure can then expel 
even more plasma from the low density region, leading to ducting instabilities. In magnetic plasma, 
e = 1 — ojp/ ( 'lo 2 — I2 2 ) for linearly polarized parallel propagating transverse waves, where 12 is the 
ion cyclotron frequency. Thus waves with w << 12 refract to high density regions, and plasma is 
attracted to regions of high wave energy density, especially when the energy and momentum of the 
waves are large. A simple expression for the ponderomotive force on an ion is derived below. 

The refraction of waves to high density regions, and the corresponding attraction of ions to 
locations of high wave energy density when the wave pressure dominates over the thermal pressure 
of the ionized component of the plasma, means that the FIP effect depends crucially on details of 
wave propagation through the chromosphere. A high wave energy density in the corona will lead to 
a strong FIP effect. Weak coronal waves and strong waves lower down may lead to an inverse FIP 
effect, as chromospheric ions are attracted downwards. We will explore precise conditions under 
which these various phenomena may occur in more detail below. 

We concentrate solely on the fractionation, referring the reader t o the r e views of Reale J2014 ) 
for the proper ties of plasma in closed coronal loops, and to Marsch ( 2006ll . Cranmer ( 2009h and 
Ofman ( 2010 1 for the physics of the acceleration of the fractionated gas into the fast and slow 
speed solar wind. In the subsections that follow, we discuss in detail the various components 
of the model. Figure [5] illustrates the basic scenario. We consider a coronal loop meeting the 
chromosphere at both footpoints. For the purposes of calculating the Alfven wave propagation, we 
ignore the loop curvature. A variety of wave processes can occur at each footpoint, some illustrated 
at one and some illustrated at the other in the figure for clarity. A steady evaporative flow with 
speed much less than the local Alfven speed is assumed to take material from the chromosphere 
to the corona. Such a flow arises from the chromospheric response to cor onal heating, as heat is 
conducted downwar ds increasing the chromospheric temperature (e.glBrav e t all 119911). Imada 


and Zweibel ( 2012h calculate evaporative upflows for a variety of coronal heating scenarios, and 
find upflows of order one to a few km s -1 in the chromosphere, similar to previous work ( Warren 


et ai. 120021: iRrav et all Il99lh . We ignore transient effects, although in a magnetic filament heated 
episodically, these could play an important role. For instance, if Alfven waves released by a coronal 
heating event reach the chromosphere before the heat conduction front, the chromospheric flow 
velocity under which the fractionation occurs will be much smaller than these estimates. In the 
opposite case, of heat conducting down before the Alfven waves arrive, no fractionation could occur. 
More exotic cases where the Alfven w aves t hemselves carry the heat that causes the evaporation, 
for example in solar flares ( Haerendel . 2009) are not considered here. 

In the absence of the ponderomotive force, we assume that the chromosphere is unfractionated, 
i.e. that sufficient turbulence exists to inhibit any gr avitatio na l settl ing of other forms of diffusion 
that might occur. A similar approach was taken bv lHansteen et aZ (1997). We suggest that the 
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reflection and refraction of acoustic waves at discontinuities associated with chromospheric shock 
waves set up the conditions necessary for a turbulence. Interactions between oppositely directed 
waves lead to fluctuations of smaller and smaller size scales, right down to microscopic dimensions 
where mixing can occur. Following on from the discussion above, we neglect diffusion processes. 
Although the thermal force is expected to develop over a similar region to the ponderomotive force 
(the region of strong density and temperature gradients in the upper chromosphere), the accelera¬ 
tion associated with the thermal force is of order 1% - 10% of the ponderomotive acceleration, and 
is negligible. The thermal acceleration is also dependent on 1/A, where A is the element atomic 
mass, unlike the ponderomotive acceleration which is mass independent, and can never produce 
an inverse fractionation as seen for example in the coronae of M dwarfs. 


6.2 The ponderomotive force 

An eyDKSsioiifoiAhejDonderoinntiveJoixejs^denved.asJollows^jsdate^fi^niAjDjDendix^in.Lam- 


ing ( 20091 1. Consider the Lagrangian for a system of n particles with mass m and electromagnetic 
waves in a box of unit volume 


r 1 / r ,2 V - '' e , r . eSE 2 — (SB 2 

L = z^ 7j m i ( v th,i + <5vi) + 2^ — (v t h,i + s Vi) • 6A + 


87T 


(1) 


where v t h,i is the thermal speed and Svi is the oscillatory speed induced by the wave of particle 
i, with mass rrii , and charge qt. Wave electric and magnetic fields are given by <5E and (SB, 
respectively, (SA is the wave vector potential, and c is the speed of light. We have chosen the 
radiation gauge where the electrostatic potential (f> = 0. In any case, this is constant in an 
electrically neutral plasma (in the absence of electrostatic waves). 

For Alfven waves, energy is partitioned according to SB 2 /8n = ]P,- mv 2 sc i /2 + SE 2 /8n. Also 
(SB • B = 0 where B is the ambient magnetic field, and Svi ■ 5A = 0, not just on average but for 
all time. We also take the time average of v t h,i • <S= 0 and v t h,i • <SA = 0 to find the Lagrangian 


i 


1 

—mv 

2 


2 

th,i 


(e - 1) SE 2 
Sn 


£ + £ 


q?SE(zi ) 2 
2 rm (fl 2 - w 2 ) ’ 


( 2 ) 


where in the second line we have written e — 1 = u 2 / (fl 2 — w 2 ) = 47mg 2 /m (fl 2 — w 2 ) and 
(e — 1) SE 2 —> Y,i4:Trq 2 5E ( Zi) 2 /m; (fl 2 — w 2 ) for linearly polarized parallel propagating transverse 
waves. The “z” Euler-Lagrange equation for particle “i” gives 


d . , q 2 d SE 2 

dt {mVth ’ iz) ~ 2m t dz (n 2 -zJ 2 ) , 


(3) 


the same as the expression derived by 

Landau et al. (19841). and agrees with earlier work fe.e 

and Parks. 19831; Li and Temerinl. 1993 

) if SE 2 = SE 2 / 2, where SE p is the peak electric field in 


wave, giving a ponderomotive force 


Fi = 


Qi 


dSE p (zi 


4 rrii (f l 2 — to 2 ) dz 


(4) 


When u> f2,, the ponderomotive acceleration is thus independent of ion mass, which is one 
crucial property relevant to obtaining an almost mass independent fractionation as observed. It 
is also hide 2 endent_onon_chax^e 1 _so_jong_j^_Jhe_^on_^sjdiax£edJand_jio^jieutraIJ ; _Litwin and 
Rosner l?1998i l give a similar expression derived from the j x B and other second order terms in the 
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Figure 5: Sche matic d i agram of model loop and wave processes, adapted from ( Laming! . 1 2012 1, 
which follows Hollweg (19841)• All footpoint wave processes may happen at both footpoints, but 
are here split between the two for clarity. Alfven waves shown as thick solid lines are assumed to 
be generated inside the coronal portion of the loop, and to bounce back and forth from the loop 
footpoints, with a probability of leaking out and being transmitted deeper into the chromosphere 
at each bounce (shown on left hand side). Approximately equal amplitudes of waves propagating 
in each direction result. Reflecting Alfven waves can also generate slow mode (here labelled as “p- 
rnode”) waves by a parametric process (shown on the right hand side as thin dashed lines). Other 
acoustic waves (“p-modes” in helioseismology parlance) propagating within the solar envelope can 
mode convert to fast mode waves upon reaching the chromospheric layer where the sound speed 
and Alfven speed are equal (approximately where the plasma /3 = 1). The upgoing fast mode waves 
(shown as thin solid lines) are refracted back downwards in the chromospheric region where the 
Alfven speed increases with height. Fast mode wave may also mode convert to Alfven waves and 
then propagate up to the loop to be reflected or transmitted, depending on the match between their 
frequency and the loop resonance. This provides an alternative means of seeding the loop with 
propagating Alfven waves. Our non-WKB wave propagation calculations in fact follow an Alfven 
wave being injected at one loop footpoint, being reflected or transmitted into the corona, and the 
successive bounces it undergoes. For waves at the loop resonance, this behaviour (without growth 
or damping included) is indistinguishable from waves generated within the loop itself, necessarily 
at the loop resonance or its harmonics. Reproduced by permission of the AAS. 
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MHD momentum equation. Away from the low frequency limit, circular polarized waves may give 
slightly different forces, with left and r ight circular polarization being different from each other as 
well as linear polarization fe.g.. lNekrasov and Fevgiil 2013 1. In nonuniform B in the low frequency 
limit, the ponderomotive force is given by 


Ft = 


4 dz 


SE P (zj) 

B ( Zi ) 2 


(5) 


which agrees with the first two terms in Equation 2.6 of lLundin and Guglielmi ( 2006i) . 

In fast mode waves, energy is now partitioned according to SB 2 /8tt + 5E t h = mv osc J 2 + 
SE 2 /8tt, where 5E t h is the wave thermal energy arising from compression and rarefaction. For the 
fast mode wave <5B • B = 0 only on average, but not instantaneously, though <5v, • S A = 0 remains 
unchanged. The thermal energy also time averages to zero, leaving the Lagrangian and pondero¬ 
motive force described by the same expressions as above. However the extra longitudinal pressure 
associated with obliquely propagating fast mode waves will reduce the eventual fractionation. See 
Sections 16.51 and 16.61 below. 

With slow mode waves, energy partitioning is similar to that for fast mode waves. Whereas 
Alfven and fast mode waves always have 5v _L B, slow modes have <5v||k. In parallel propagation 
the electric field, and hence a ponderomotive force capable of ion-neutral separation, are absent. 
Transverse electric field may arise with a_slow m ode wave in paralld propagation along an isolated 
flux tube (a “sausage” mode) (e.g., Mikhalvaev and Solov’ev! . 2005!). We do not consider such 
possibilities further. 


6.3 Chromospheric model 

The ponderomotive force depends on the gradient of the wave transverse electric energy, d8Ej_/dz. 
Since SE± = 5v±B/c, such a gradient will develop wher e Sv \ varies in response to varying chromo¬ 
spheric density. Previous (mainly analytic; e.g. iHollwed . 119841 works on wave propagation through 
the chromosphere have approximated the vertical density structure as exponential, with a scale 
height of order 200 km. While adequate for studying wave transmission into the corona, more 
detail is required for the FIP effect. The chromosphere is heated, presumably by means similar to 
those thought to heat the corona, and cooled in the region where most FIP fractionation occurs 
principally by radiation in H Lyman a. In this region where H starts to become ionized, its radia¬ 
tive cooling becomes increasingly inhibited and the temperature must rise. Correspondingly, the 
density drops, and the density gradient so produced is steeper than the typical hydrostatic scale 
height elsewhere. Consequently, the region of the chromosphere where the FIP fractionation has 
long been assumed to occur, i.e. the region where H becomes ionized, is highly likely to be a region 
with a strong density gradient due to the physics of radiative cooling, and hence a location where 
the ponderomotive force will be stron gest. F igure [5] sh ows this structure in the left panel, taken 
from the empiri c al mo del (C7) of Avrett and Loeser ( 2008h . an update to the older VALC model 


( Vernazza et al . Il981 ). 


The right panel of Figure [6] shows the model magnetic field. Altitude here is given by the y-axis. 
The x-axis gives lateral expansion on the same scale. The dashed lines show contours of the Alfven 
speed. The magnetic field emerges from the photosphere in tight fibrils. At the atmospheric layer, 
typically in the low chromosphere, where gas pressure and magnetic pressure are equal, the field 
begins to expand and fill the whole volume. In Figure [6] (right panel) this corresponds to y = 650 
km approximately. In the real Sun, this occurs at altitudes typically between 400 and 800 km, 
depending on the magnetic field strength. The total magnetic fiel d expa nsion from y = 600 km 
upwards in this model is a factor of 5, similar to that suggested bv iGar ^ (l200lh. 


Observations of hard X-ray bremstrahlung during flares with RHESSI jLin et al . 2004 1 cor- 


V/UOOl V CAiUXVXXXO VXX XXCAX VX J. A X <AJ L»X OXXXO LX Uiii LIXXp, VX LAX 1XXL, XidX OO VVXUXX 1L11UUU1 \ XJXXX U IXX . . AjVJVXT: J V_,VXX 

roborate many of these features ( Kontar et al . 20081 : (Saint- Hilaire et al . 2010l f. An average scale 
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Figure 6 ^^efc^m£iricaHiiode^fthe_solar_chroinos 2 liere - ^lotted_^oin - data i _giveiiJr^Avrett and 
Loeser lf2008l ). The density is the solid black line, to be read on the left axis, while the temperature 
is the dashed gray line, to be read on the right axis. The ordinate (“distance along the loop”) 
is the altitude above the photosphere. The steep density gradient at 2150 km altitude is where 
strong ponderomotive f orce c an develo p. Right: For ce-free magnetic field structure calculated in 


Rakowski and Laming ( 2012 ) following Athav ( 198ll ). Altitude is given by the y-axis. The x-axis 


give lateral expansion. The dashed lines show contours of the Alfven speed. In this plot the plasma 
/? = 1 layer is at an altitude of approximately 650 km in the Sun. Reproduced by permission of 
the A AS. 


density scale height of 130- 140 km is derived. Kontar et al. 120081 ) infer chromospheric flux tube 
expanding in radius by a factor of 3 at an altitude between 900 km and 1200 km. This would 
imply magnetic fie ld de creasi ng by a factor of 9, albeit at a higher height than suggested above. 

Avrett and Loeser ( 2008 ) also give an empirical electron density coming from the ionization 
balance for H, from which we calculate the ionization balance for all other elements of interest. 
Collisional processes (io nizati on, radiativ e and dielectronic recombination) are all taken from the 
tabulations provided bv iMazzotta et ~ai. ( 19981) . Subsequent ref inem ents to the ionization and 
dielectronic recombination rates as implemented bv lBrvans et all (120091 ) are not crucial here, since 
our concern is only with neutrals and singly ionize d spec ies, but the suppression of dielectronic 
recombination by finite density e ffects (Nikolic et all. 12013 ) is important and is now included as an 
update to previous c alculations (Laming . 20091.12012 ). We also include charge exchange io nization 
using rates given by iKingdon and Ferland ( 19961 ). updated and amended as detailed in Ibamingl 
( 2012 ?). This treats the capture by a free proton of an electron from an initially neutral atom. 
It is most important for O, and helps lock the O ioni zation balance to t h at o f H, but can be of 
significance for other low FIP elements as well. In the lAvrett and Loeser (2008 ) model, H begins 
to become ioniz ed at he igh ts betwee n 1000 1500 km, reaching about 50% ionization at 2000 km. 
For comparison, ISaint-Hilaire et al. ( 2010l ) infer a unit step change in ionization level at a height 
of 1.3 ±0.2 Mm. 

We evaluate photoionization rates using incident spectra b ased onlVernazza and Reeves! ( 1978h 
with the extensions and modifications outlined in Laming! ( 2004a ). In most cases, the “quiet 
region” spectrum is used. An additional contribution from trapped Lyman a in regions where 
this line is optically thick is also included, with approximate flux f nn dlAC ex / {A + rCde-ex) 
where J dl is the column density of H atoms, r = a f ?ih dl is the opacity in terms of the 
absorption cross section a at line centre, A is the upper level decay rate, and C ex and Cde-ex 
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are the collisional excitation and de-e xcitatio n rate s for the upper level. We take photoionization 
cross sections from the compilation of Verner et al. ( 19961 ). Our ionization balance so computed is 
necessarily a steady state model, but is based on the empirical electron density found bv Avrett 
and Loeser ( 2008t h which is highe_r_ than that which we would compute based on their inferred 
chromospheric temperature. Carlsson and Stein! (1 2002 1 study the effects of dynamic ionization 
of H, and find that the passage of shock waves through the chromosphere does indeed increase 
the electron density over that coming from simple photoionization-recombination equilibrium for 
H. Individual photoionization and recombination rates however are too slow for the ionization 
balance to respond to individual shock waves. The elevated electron density instead represents a 
time av eraged response to the dynamics, and is itself ouasi-steadv state. Wedemever-Bohm and 
Carlsson (12011 1 consider the Ca + to Ca 2+ ionization balance. It is more variable than that between 
H and H + , but is of less concern to us, since the ponderomotive force experienced by an ion is 
independent of ion charge, so long as the wave frequency is much lower than the ion gyrofrequency. 

In isolation, the action of the ponderomotive force would produce a local abundance variation 
close to the region where it is strongest, around chromospheric altitudes of 2150 km. A region 
where low FIPs are depleted by the upwards acceleration would sit below this altitude, and above 
it a region of FIP enhancement. A coronal abundance anomaly would not result without either a 
flow through the chromosphere, or diffusion (in a static case). The typical timescales associated 
with such processes would be Hjj/v ~ 2 x 10 7 /10 3 ~ 2 x 10 4 s where Hp is the density scale height 
and v the flow velocity lower down in the chromosphere, or H^/D ~ 4 x 10 14 /6 x 10 9 ~ 6 x 10 4 
s where D is the diffusion coefficient for low FIP ions in a neutral H atmosphere . Both timescales 
evaluate to several hours to days, similar to that observed ( Widing and Feldmanl . 2001 1. as do 
similar timescales evaluated for the coronal loop. The chromospheric response may be altered by 
the mixing induced bv turbulence associated with shocks in the lower chromosphere (e.g.. Reardon 


et al., 20081 ). and so the coronal diffusion is likely to be the controlling timescale. We do not 
consider such processes further, restricting ourselves to a steady state chromosphere and upward 
flow speed, and emphasizing again that collisional coupling in the chromosphere is crucial to the 
coronal abundance anomaly, since it is the source for the extra low FIP ions. Once collisionless 
(in the solar wind) no further bulk fractionation can occur. A local FIP enhancement must be 
accompanied by another local depletion, since the collisional coupling to the lower solar atmosphere 
has been lost. 


6.4 The Alfven wave transport equations 

Alfven wave propagation is often treated in the Wentzel-Kramers-Brillouin (WKB) approxima¬ 
tion, where it is assumed that the wavelength is much shorter than the typical length scale of 
inhomogeneities, and hence that the effects of reflection and refraction can be neglected. This is 
emphatically not the case in the solar chromosphere. Since the ponderomotive force depends on 
the spatial variation of the Alfven wave electric field in the inhomogeneous plasma of the chro¬ 
mosphere, a direct integration of the transport equations is required. We evaluate this spatial 
variation paraphrasing and updating the treatment in Appendix B of lLamin 5 (120091 ). We start 
from the linearized MHD force and induction equations, 


d5v 

p-gjT + P ( u ' V) 6v + (6v ■ V) pn = 


(V x <SB) x B (B • V) <5B - (V<5B) • B 


47T 


4-7T 


( 6 ) 


and 


dSB 

~df 


= V x (<Sv x B) + V x (u x (5B) 

= —BV • 6v — (<5v • V) B + (B • V) <5v — <5BV • u — (u • V) <5B + (<5B • V) u, (7) 
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where u and B are the unperturbed velocity and magnetic field, <5v and <5B are the perturbations, 
and p is the density. Equation © is rewritten using (V5B) ■ B = V (B ■ (5B) — (VB) • <5B to yield 


dSv 

—— + (u • V)6v = V A • 
at 



£B V A -Vp (VB) • £B _ (<W ■ V) (pu) 
y/Airp 2 p Anp p 


( 8 ) 


where V a = B / s/Anp is the Alfven velocity, and we have assumed B ■ 5B = 0 for Alfven waves 
and parallel propagating fast modes. 

Specializing to plane Alfven waves, we write (VB) • <JB = (dB x /dx)S B = — (dB z /dz) SB/ 2 
since V ■ B = 0 (assuming dB x /dx = dB v /dy in cylindrical symmetry, with <5B = SBx where x 
is a unit vector along the x-axis). Similarly, (<5v • V) (pu) = Sv x d (pu x ) /dx = —5v x d (pu z ) /dz/2 
since V • (pu) = 0, and using d ( pu z /B z ) /dz = 0 gives 


dSv 

— + (u ■ V)5v = 



<5B V A <SB Va , u 
yJAnp 2 H d \fAnp 2 H b + 2H B 


( 9 ) 


Here I/Hb = dlnB z /dz, 1 /Hjj = dlnp/dz , and below A/Ha = 9In V A /dz. 

Similar manipulations using V ■ u = — u • Vp/p = —u/Hd and assuming V • 5v = 0 give the 
induction equation in the form 


d ( SB \ . m 

\~*= +(u-V) 


SB 


= (V A • V) 5v + 


5B 


<5v' 


V A 


SB 


dt \y/Anp J ' v ' yJAitp v ^ ’ ' y/Anp 2 Hd ' 2 Hb \JAnp 2 Hb 

Taking Equation © plus or minus Equation (1101) and rearranging gives the final result, 


^ + (uTV A )^- = (u±V A ) 


AH d 2 Ha 


( 10 ) 


( 11 ) 


where I± = Sv ± SB/y/Anp, representing waves propagating in the =F z-directions. We have 
focussed on plane polarized Alfven waves. Clearly circularly polarized Alfven waves will obey the 
same transport equations, as do torsional waves in cylindrical coordinates. These different wave 
polarizations may still produce subtly different fractionations, due to their different couplings to 
other wave modes. The case of parametric generation of slow mode waves is discussed below. 

Such transport equations h a ve be en presented in several different forms by various authors. 
ICranmer and van Ballegooiienl ( 2005l l review these in their Appendix B and demonstrate their 
equivalence. We have assumed a cylindrically symmetrical magnetic field above. One could easily 
generalize this treatment to incorporate Alfven wave transport on arbitrary (i.e., observed or 
extrapolated) magnetic field lines for a more realistic description of observations. We have also 
restricted treatment to field aligned wave propagation, rendering the waves incompressible. For 
plane and circularly polarized waves this amounts to assuming that the wave follows the field line 
on which it propagates. For torsional waves on a flux tube, it requires that the flux tube magnetic 
field have no Bp component. In the case that this is not so, the torsional wave is inevitably of 
mixed Alfven/fast mode polarization and has B • dB ^ 0. In this case Equations 15 are augmented 
by extra terms —V (B • <5B) /47rp =F VaV • Sv on the right hand side. 

In integrating Equations CD. we put dl±/dt = iu)I± to derive four coupled equations for the 
real and imaginary parts of I±. Equations (HID are integrated from a starting point in the left 
hand side chromosphere where Alfven waves leak down into the chromosphere, back through the 
corona to the right hand side where waves are fed up from below. In this way the reflection and 
transmission of Alfven waves at the loop footpoints and elsewhere is naturally reconstructed. The 
velocity and magnetic field perturbations are calculated from 


Sv = 


/_ 
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The wave energy density and positive and negative going energy fluxes are 

— - - (I 2 + I 2 ) 

8n 


and the wave peak electric field appearing in equation 4 is 


u = 

pSv 2 

2 + ' 

F+ = 

Pt2 V 

4 T + Va 

= 

P t'2 t/ 
\ T ~ Va ' 


( 12 ) 


(13) 


(14) 


6.5 Fractionation 


We recap and expand the treatment given in Lamind (2004a). Our model starts with a fully mixed 
chromosphere, upon which pondermotive forces due to Alfven wave reflection and transmission 
act to provide the fractionation. The low solar chromosphere is of much higher density than the 
upper layers where the FIP fractionation will occur in our models. Consequently the lower bound¬ 
ary condition of completely mixed photospheric composition material gives an essentially infinite 
particle “reservoir” to supply the extra frac tionated element s. To calculate the fractionation, we 
follow in part the approach and notation of Schwadron et al. ( 19991) . Consider first the motion of 


ions and neutrals of element s in a background flow of protons and hydrogen with speed u. We 
neglect the ambipolar force, which is generally much less than gravity. We also neglect an inertial 
term d/dz (p s u 2 / 2) since the flow speed will turn out to be much lower than the particle thermal 
speeds. Then the momentum equations for ions and neutrals are 


dP S j 

dz 

dPsn 

dz 


— Psid Psi^si i^si ^) 

= PsnP PSnVsn (^sn ^) 5 


(15) 

(16) 


where P S i and P sn are the partial pressures of ions and neutrals of element s, p S i and p sn are 
the corresponding densities, iz si and v sn the collision rates with ambient gas (assumed hydrogen 
and protons), u s , and u sn the flow speeds, and u the hydrogen flow speed imposed on the loop. 
Detailed expressions for iz si and v sn are given in iLaminej ( 2004a ). The partial pressures are given 
by P s ,i,n = (/cB^/ms + ?)( urb + rj, ave ) p Sl i, n /2, where the various terms on the right hand side 
represent the particle thermal speed, the microturbulent velocity coming from the chromospheric 
model, and the particle motion in longitudinal (i.e., acoustic) waves, either induced by the Alfven 
waves themselves, or propagating up from the photosphere. These are discussed in more detail in 
the next subsection. In true collisionless plasma, neutrals would not respond to the wave. However 
the solar chr omosphere is sufficiently collisional that neutrals move with the ions in the wave 
motion (e.g., Vranies et all [20081 ) for wave frequencies well below the charge exchange rate that 
couples neutrals and ions, and so neutrals require the same form for their partial pressure as the 
ions. We also neglected inertial terms above. Inclusion of such terms would lead to 


— (&B T/m s + fturb T 'Wave A u s ) Psi,n /2 — V s Psi,n /2 


(17) 


assuming u S i — u sn -C u s ~ u. 
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The momentum equations are combined by forming v sn dP S i/dz + v S idP S n/dz to yield 


dP s , > , di s p s v 2 v si - v sn 

~K~ = ~Ps9 ~ VeSPs ( U s -u) + j -- 7 - 

2 (1 “1“ Ss^S' 


( 18 ) 


with i/ e ff = v S iV S n/ (JisVsn + (1 — £, s ) v s i) and £ s being the ionization fraction of element s. For 

u s = u- g/v eS (1 - n/m s ) + d^/dz (i> 2 /2) (v si - u sn ) /v ai /v sn , (19) 


where p is the mean molecular weight, all elements are lifted by the background flow to the same 

scale height given by k^T/pg. This obviously requires uv e s > g , assuming the term in d^ s /dz 

negligible. Generally u — g/v e s (1 — £t/ m s) ~ 10 3 cm S_1 or more, while ~ 10~ 8 cm -1 and 
2 

( v si — v sn ) /v si /v sn ~ 10° cmV 1 . For uv s <C g, we get gravitationally stratified solutions 
with p s oc exp (— m s gz/k-aT ). With g = 2.74 x 10 4 cm s -2 , and v s ~ 10 2 — 10 3 the former case is 
vali d for the Sun for flow speeds in the chromosphere greater than a relatively modest 10 100 cm 
Laming ( 2012 1 considers the upward flow speed due to chromospheric evaporation, and finds 


speeds of at least 10 3 cm s _1 . 

A gravitational settling timescale may be estimated from equations [15] and [T6] We write 
dP s /dz ~ Cgdps/dz ~ (cg/u se tt) dps/dt where the sound speed is cs and the settling speed 
u se tt ~ g/i'eff- Integrating with respect to t gives t ~ c 2 s v eff/g 2 ~ 10 4 {jih/ 10 10 ) /A s, where A is 
the atomic mass of element s, assuming that v e ff is dominated by collisions with neutral H. The 
settling time increases with the density as one moves deeper into the chromosphere, and with the 
presence of turbulence with amplitude added in quadrature to cs, and can easily be on the order 
of days or weeks or more. Higher in the chromosphere, v e ff will be dominated by collisions with 
protons, and a longer timescale would result as the chromospheric density approaches 10 10 cm -3 . 

The solar chromosphere is doubtless a more dynamic environment than represented by Equa¬ 
tions (USD CEO). For our purposes the net result of this dynamic behavior is assumed to completely 
mix up the plasma to give uniform elemental composition with height, which is obtained in our 
model with the above choice for u s . Other choices may be possible which would provide chemical 
fractionation in the unperturbed chromosphere, and one could choose u s to provide the required 
FIP effect. However the physics behind such a specification for u s in most cases remains obscure, 
and is probably unrealistic, leading to an unsatisfactory explanation for the FIP effect. 

We now include a ponderomotive force, p S ia on the ions in the momentum equations; 


dP S i 

dz 

dPsn 

dz 


— Psig PsiVsi {Vsi u) T PsiOj 

= Psn.g Psn^sn (u S n ^) • 


( 20 ) 

( 21 ) 


Again writing P s = p s c|, taking u s as specified above and integrating with respect to z we find 


Ps {Zu) 
Ps (zi) 


= exp 


h; 




dz 


( 22 ) 


where the constant of integration has been chosen to keep p s (z u ) = p s ( zi ) when a = 0. A quanti¬ 
tative assessment of coronal element abundances anomalies requires an evaluation of Equation (1221) 
with a realistic model chromosphere in the region of Alfven wave reflection, with a coming from 
Equation (O. 

As will be seen below, most the the FIP fractionation develops over a small range in z, where the 
ponderomotive acceleration a is strong. The timescale to establish this local abundance anomaly is 
short, basically (z u — zi) /u ~ (l0 6 — 10 7 ) /10 4 ~ 10 3 s. By itself, the ponderomotive force would 
only produce a deficit of low FIP ions just below this region, and a surplus just above. To produce 
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the fractionation seen in coronal loops and the solar wind, thermal transport or diffusion processes 
must continually supply further low FIP ions from below to erase the deficit, and communicate the 
resulting fractionation above to still higher levels of the solar atmosphere. Thus while the fraction¬ 
ation can be produced locally rather quickly, and is therefore immune from processes that would 
otherwise disrupt fractionation produced solely by diffusive processes, the timescale to change the 
composition of a coronal loop is still on the order of days, as observed ( Widing and Feldmanl . 2001 1 
due to the necessity of thermal transport of the fractionation. 


6.6 Compressional chromospheric waves 

6.6.1 Introduction 

We have introduced an extra longitudinal pressure associated with the Alfven waves proportional to 
v^ ave (see Equation (fT71) ), which has the effect of causing some saturation of the FIP fractionation. 
Here we give some physical motivations for this extra term. It has three possible sources. The 
first is the inevitable generation of slow mode waves by the Alfven driver. Physically, the periodic 
variation in magnetic pressure of the Alfven wave drives longitudinal compressional waves. These 
generated acoustic waves can act back on the Alfven driver, as the compressional wave introduces 
a periodic variation in the Alfven speed, which generates new Alfven waves. The second will 
be obliquely propagating fast mode waves, which are necessarily compressive. Torsional Alfven 
wave on a twisted magnetic flux tube are inevitably mixed with the kink mode, and hence are 
compressive. The third relates to acoustic waves propagating up to the chromosphere from the 
convection zone. 


6.6.2 Parametric generation by Alfven waves 

Following iLamina ( 2012 1 we illustrate the generation of slow mode or acoustic waves by the pon- 
deromotive force associated with plane Alfven waves with a simple ID calculation. The linearized 
momentum equation keeping terms to all orders in perturbed quantities is (all symbols have their 
usual meanings as defined above), 


, e , dSv z , . , dSv z , . 

(P + Sp) + (p + dp) Sv *-qT = (P + 5p ) 


d SB 2 dSP 

dz 87r (p + Sp) dz 


-gdp, 


(23) 


where 

6p = -pV • dr - Sr z |f = -pSr (ik s + 

SP = -yPV • <5r - Sr z = -PSr (iks'y + , (24) 


for Eulerian displacement Sr oc exp* (wgt + k$z), where Hp = P/ ( dP/dz ) and Hp> = p/ ( dp/dz) 
(signed) pressure and density scale heights, respectively. The first term on the right hand side 
of Equation (l?51) represents the instantaneous ponderomotive force. In cases where the WKB 
approximation applies, SB oc p 1 / 4 , and this expression is equivalent to the more usual form 
—d (5B 2 /8n) /dz. Substituting for Sp and SP from Equations (l24l) . keeping terms as high as 
second order in small quantities, Equation (l23l) becomes 


■ W ‘ s ' s 2 i ( 2,722 

- i-j-j—8v z -I- I -u} s + k s c s - i 


r 2 

c s 


lH % - ‘w; - ikss - t) Sv • ■ 


” s !it =0 - (25) 


This is considerably simplified in isothermal conditions, 7=1, Hp = Hp> = —c 2 s /g , so that 

. iOg c n / 2 i22 -7 \r ■ S) SB 

- l Jj^ Sv z + + k s c s + iksg) Sv z - luJ s-gZ-^ = 0. 


(26) 
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We put Ssks = —g/2cg, and = 2SftAu/ra, u = 2ua/ti, where kA and ua are the 

wavevector and angular frequency of the Alfven wave with n = 1,2,3,.. (anticipating the result 
below). We find 


5v z - 5v z iH D u s 1 - 


n 


IID 


d SB 2 

dz 8np 


= 0 


with solution 


6v z = 


—i 


'Sv 2 a 


+ H 


d m s 


1 - 


VJ 


- Hdws 1 — 


Vi 


(27) 


(28) 


where we have put (5B 2 /&Trp) = (di? 2 /47rp) /4 H D = 8v 2 A /^H D using the WKB result for 
Alfven waves in a density gradient, and assuming 1 /H D When c| ~ V\ or H D 0, 

5v z ~ —i8vA/2. In the opposite limit Sv z ~ — i5v 2 A /AHduj s (l — Cs/^a) ■ I* 1 these two cases 
us = u>a or u>s = 2 u A , respectively. In fact, acoustic waves can be generated with us = 2 uja/ti, 
with higher n becoming mo re intense as t he nonlinearity increases ( Landau and Lifshitzj . Il976lf . 
This is shown explicitly by Laming (201 2h. wh o extends the treatment above to include higher 
powers of Sv z . IVasheghani Farahani et all ( 201 il l treat the case of slow mode wave generation by a 
torsional Alfven wave. This is subtly different to the case of a plane Alfven wave considered here, 
and the FIP fractionation resulting from such a model will be investigated in future work. 

When the Alfven wave becomes evanescent in a sufficiently steep density gradient, kA —> in A 
and 


Sv z = 


—i 


l5v 2 A 


II 


D^S 


1 + 


,2 2 

X A C S 


U“, 


— Hdus 


2 2 \ 

1 + ^) 
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resulting in reduced 5v z compared to the case with propagating Alfven waves. This will be apparent 
in results to be presented in Section [T] 


6.6.3 Fast-mode waves from above 


The treatment of Vasheg hani Farahan i et aL ( 20111 ) also considers a second source of longitudinal 
pressure. Torsional waves on a flux tube with magnetic twist are no longer pure Alfven waves, but 
have an admixture of fast mode or kink polarization. Since the magnetic perturbation (in the <t> 
direction in cylindrical polar coordinates) now has a component along the background magnetic 
field (the twist gives it a (f> component also), the wave becomes com pressiv e. This is also true 
more generally of any fast mode wave in oblique propagation, and in lLamin 5 (l2009t) it was shown 
that such longitudinal pressure would lead to a natural saturation of the FIP effect at high wave 
amplitudes. In the case where transverse and longitudinal pressures are equal, the FIP effect 
saturates at the commonly observed value of about 4. However such compressive waves are also 
subject to increased damping, and are therefore less likely to develop to large amplitudes. The 
existence or otherwise of twisted magnetic loops in the corona is also controversial. Accordingly, 
we do not consider such possibilities further in this review. 


6.6.4 P-modes from below 

The third contribution to the longitudinal pressure in the chromosphere comes from the continued 
propagation upwards of slow modes waves or shocks developing from turbulence lower down in 
the atmospher e. We include a longitudinal pressure with amplitude 6.25 km s _1 motivated by the 
simulations in Heggland et al. (|201ll ). The velocity amplitude of these waves would be expected 
to grow with altitude, due to WKB effects, although some reflection back downwards may also 
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occur, and should damp mainly due to the effect of radiation. The principal co olant in t he p hoto- 
spher e is H _ , with cooling timescale r ~ (l.2 x 10 1:L /n e ) (5000 K/T) 2 minutes (Ayres and Rabin. 


1996). Cooling of a similar magnitude will be present in the low chromosphere, although different 


processes may also contribute. Although the plasma density varies, throughout large sections of 
the chromospheric model, the electron density and temperature are roughly constant and of order 
10 11 cm~ 3 and 5000 K, respectively, leading to a constant cooling time of order 1 minute. Thus for 
3 minute or 5 minute waves, t he product uit ~ 1 , and the wave decay length would be comparable 
to the wavelength jMrhajug_ajrdJWeibel-Mihalasl . Il999h . and similar to the WKB growth. In fact, 
the simulations of lHeggland et all ( 2011 1 show approximately constant slow mod e wave amplitude 
with height for a variety of models, as do the observations of lBeck et al. ( 20131 . albeit at a lower 
wave amplitude), so we accordingly adopt a constant slow mode wave amplitude in our models. 

While slow mode waves reduce the fractionation by increasing v1 in the denominator of the 
integrand in Equation (1221) . when these steepen into shock waves, we expect fractionation to cease 
altogether due to the extra turbulence. Accordingly we also neglect fractionation at chromospheric 
altitudes where the slow mode wave amplitudes added in quadrature come to a value greater than 
the local sound speed. The rationale for this is that in such conditions, a shock will form, and the 
mixing produced by turbulence generated in such conditions would preclude further fractionation. 
Reflection and refraction of sound waves at the discontinuity will produce oppositely directed 
packets of sound waves which can initiate a cascade to smaller, and ultimately microscopic, spatial 
scales, assumed to induce microscopic mixing. This has the effect in most models of limiting FIP 
fractionation to the upper regions of the chromosphere. 
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7 Results and Interpretation 


7.1 Closed loop 

As a first illustration of the action of the ponderonrotive force in causing FIP fractionation, we 
consider the case of a closed magnetic loop. The model loop is taken to be 100 000 km long, 
and plasma in the loop has a density scale height of 125 000 km, extrapolated upwards from 
the chromospheric model. The coronal temperature is constant at the maximum temperature 
of the chromospheric temperature (~ 1.6 x 10 6 K) but this value is not significant since we are 
uninterested in the ionization balance in this region and the Alfven wave propagation only depends 
on the density. The coronal magnetic field is 20 G, increasing by a fac tor o f 5 i n the p hotosphere, 
as in Figure [6] (right panel). The chromospheric model is taken from Avrett and Loesen ( 2008 1. 
and the plasma /? = 1 layer is at 615 km altitude above the photosphere. We simply treat the 
parallel propagation of undamped shear Alfven waves. With reference to Figure [SJ the solution 
of Equations m is started with a downgoing wave at the left-hand chromosphere, and then 
integrated back to the chromosphere at the other end of the loop. Choosing a wave angular 
frequency of 0.07 s -1 means the wave is on resonance with the coronal loop, and in our calculation, 
98% of the upgoing wave energy from the right-hand chromosphere is transmitted to the corona. 
Hence the chromospheric wave patterns in both chromospheres are very similar, with the roles of 
upgoing and downgoing waves interchanged (with 100% transmission, they would be identical). 
The ponderomotive force is insensitive to the direction of wave motion, and so is the same at each 
footpoint. 

Figure 0 shows the solution for the coronal portion of the loop. The top panel shows <5v 
and 5T5/y/4:irp (Equations Ijl2[l ). real components in black, imaginary components in grey. The 
middle panel shows the wave energy fluxes, left and right going, and their difference, divided by 
the magnetic flux density. Normalized thus, this difference should be constant in the absence of 
wave growth or damping, and provides a check on energy conservation in the integration. The 
bottom panel shows the ponderomotive acceleration, a, calculated from Equations © and m- 
Throughout the corona (2500 km < z < 97500 km), a <C 2.7 x 10 4 cm s 2 , the solar gravitational 
acceleration. In the chromospheric footpoints, a is much larger. 

Figure [5] illustrates important features of the solution in the chromospheric portion of the 
loop. The top and bottom panels on the left hand side show the same variables as the top and 
bottom panels in Figure [7] (the wave energy fluxes are not shown). The ponderomotive acceleration 
shows a large “spike” at an altitude of 2150 km, coinciding with the steep chromospheric density 
gradient shown in Figure ED (left panel). Also shown on the bottom left panel is the amplitude of 
slow mode waves arising from parametric generation by the Alfven driver, shown as a dotted line. 
This is added in quadrature to the amplitude of slow mode waves propagating upwards from the 
photospheric, assumed here to be a constant 6.25 km s~ 4 , following the discussion in Section I77T7T1 
T his value is also consi stent with the macroscopic velocity field inferred in the solar chromosphere 
by Vievtes et al. ( 20051 . see their Figure 3ab). The parametric slow mode amplitude goes through 
a sh arp minimum wher e the c hromospheric density gradient is steepest at about 2150 km altitude. 
iRakowski and Lamind ( 2012l l show that this occurs where the Alfven wave becomes evanescent. 
The panels on the right hand side show the chromospheric ionization balance and FIP fractionation 
for selected elements. The bottom panel shows FIP fractionation for the abundance ratios S/O, 
He/O, Mg/O, Fe/O and Si/O, with the linestyles indicated as black lines, to be read on the left 
axis. One can see Fe/O and Mg/O fractionating to a value close to 4, Si/O to 2.6. S/O to 1.6 and 
He/O to 0.4. The ionization balance is shown as grey lines with the same linestyles for Fe, Mg, 
Si, S, C, and He. Also shown are charge states fractions for H and O as labelled. The ionization 
balances for Fe, Mg, Si, C, and S are hard to distinguish on this plot, so the top right panel shows 
these in more detail, with the charge fraction scale running from 0.99 to 1.00. 
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Figure 7: Solution for the coronal portion of the loop. The top panel shows <5v and (5B/y / 47rp 
(Equations (1121) '). real components in black, imaginary components in grey. The middle panel 
shows the wave energy fluxes, left and right going, and their difference, divided by the magnetic 
flux density. Normalized thus, this difference should be constant in the absence of wave growth or 
damping, and provide a check on energy conservation in the integration. The bottom panel shows 
the ponderomotive acceleration, a. 
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Figure 8: The top and bottom panels on the left hand side show the same variables as the top 
and bottom panels in Figure 0 (the wave energy fluxes are not shown). The panels on the right 
hand side show the chromospheric ionization balance and FIP fractionation for selected elements. 
The bottom panel shows FIP fractionation for the abundance ratios S/O, He/O, Mg/O, Fe/O and 
Si/O, with the linestyles indicated as black lines, to be read on the left axis. One can see Fe/O and 
Mg/O fractionating to a value close to 4, Si/O to 2.6. S/O to 1.6 and He/O to 0.4. The ionization 
balance is shown as grey lines with the s ame linesty les f or Fe, Mg, Si, S, and He. He stays neutral 
the longest, and matches calculations in lGolding et all ( 2014 ). Also shown are H and O as solid 
lines, labelled respectively. The O ionization balance follows that for H due to the strong coupling 
by charge exchange. The ionization balances for Fe, Mg, Si, and S are hard to distinguish on this 
plot, so the top right panel shows these in more detail, with the charge fraction scale running from 
0.99 to 1.00. 
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Fe and Mg have similar charge state fractions for all altitudes of relevance, and hence fractionate 
to about the same degree. Si, although low FIP and 99.8% ionized in the region 1500 km to 
2500 km, fractionates noticeably less. Around 1000 km altitude, Si is more highly ionized due 
to the increased flux in Ly a , which is now trapped by the increased opacity in the line. In the 
model presented here, no fractionation occurs at this height. S is also has a high charge state 
fraction, at 99.4%, but its ionization potential is higher than the energy of Ly a, and so it is not 
affected by the trapped radiation. This dependence on the ionization fraction comes from t he fo rm 
of v e f f = v S iV sn l {i, s v S n ± (1 — £ s ) v s i) in the integrand of equation 1221 where v 8 i » u sn , I Lamind . 
l2004af) . which is the case when the background gas is mainly ionized H. Only when ( s ~ 1 to a 
high degree of precision does the second term in the denominator become negligible compared to 
the first, allowing strong fractionation to occur. Lower down in the chromosphere, where neutral 
H dominates, i/ sn ~ v S i , and this precise dependence on element charge state is lost. 


Table 3: FIP Fractionations in Closed Magnetic Field (see text for details). 


ratio 


models 




observations 





75 

95 

110 

a 

b 

c 

d 

e 

f 

g 



( km s _1 ’ 

) 








He/O 

0.56 

0.40 

0.22 

0.68-0.60 

0.54-0.98 





0.29 

C/O 

1.14 

1.15 

1.03 

1.36 1.41 

1.23-1.55 





0.76 

N/O 

0.85 

0.73 

0.55 

0.72 1.32 

0.59-0.68 





1.02 

Ne/O 

0.71 

0.55 

0.36 

0.38-0.75 

0.79- 1.15 





0.75 

Na/O 

2.57 

3.95 

5.83 


2.51-3.31 

1.8+^ 

7.8+g 3 


3.43 ± 1.7 


Mg/O 

2.37 

3.52 

4.95 

2.58-2.61 

1.95-3.55 

2.7T0.3 

2.8+ 3 



2.71 

Al/O 

2.20 

3.16 

4.28 


2.29-3.02 

5.6+^ 

3-6+2 




Si/O 

1.92 

2.60 

3.26 

2.49-3.11 

2.14-3.26 


4 -9+s 


2.5±1.0 

2.37 

P/O 

1.52 

1.83 

2.01 



c +11 
°-3.4 




S/O 

1.37 

1.56 

1.61 

1.62 1.92 

1.20-2.09 

2.1±0.2 

2.2T0.2 

1.7 ±0.3 

L°° +32 

1.00 

Cl/O 

1.04 

0.97 

0.84 






1.8±1.8 


Ar/O 

0.93 

0.83 

0.65 





1.1 ±0.1 

1.00±0.15 

0.73 

K/O 

2.71 

4.29 

6.58 



4 7 + 7 - 0 
' - 2.8 

1 &+ 0,4 
- 0.6 

3.5 ±0.9 

c c +9 

0,0 -3.7 


Ca/O 

2.71 

4.29 

6.58 


2.09-3.88 

2.7T0.25 

o c +4.3 
00 -1.9 


3.0-9.7 

2.86 

Ti/O 

2.70 

4.28 

6.57 







Cr/O 

2.66 

4.18 

6.36 


2.40-3.47 






Fe/O 

2.66 

4.18 

6.35 

2.28-2.90 

1.95-3.55 


7 o +1,4 
' -1.2 



2.27 

Ni/O 

2.39 

3.57 

5.11 








Kr/O 

0.99 

0.89 

0.73 









Table [3] gives model fractionations for three different assumed coronal Alfven wave amplitudes, 
75, 95, and 110 km s+ which are initiated with amplitudes 0.6, 0.75, and 0.9 km s _1 low in 
chromosphere A. Figures [7] and [8] correspond to the middle entry, 9 5 km s -1 . The model results are 


compared with obs e rvation s taken from the lit erature; (a) Zurbuchen_et al J 20021 ) . given relative 


to O, (b) Bochsle~rl ( 2007al ). relative to O, (c) iGiammanco et aH ~( 20071 . 20081). relative to H, (d ) 


Brvans et al. (2003), given relative to the me an of O, Ne and Ar and, (e) Phillips et al\ d2003h 
relativ e to H, (f) Ar, Svlwester et al. (12010a); K, Svlwester et al. ( 2010bl): Cl. ISvlwester et all 


( 2011 1; S, Svlwester et al. ( 2012 1: Si, Svlwester et al. ( 20131 ); an d Na Phillips et al. ( 2010 1. all 


relative to H and given r elative to the p hotospheric abundance of Asplund et al\ (120091 1 , and (g) 
SEP m easu rements from Reames l (1 201 ill . given relative to O. Zurbuchen et al. ( 2002 1. Bochsler 
(2007a) and Giammanco et al. ( 200*^20081 1 give abundances derived from in situ measurements in 
the slow solar win d. Th ey genera lly agree best with model with coronal Alfven wave amplitudes of 
75-95 km s _1 . Brvans et al. ( 20091 ) give abundances derived from a deep SOHO/SUMER (Wilhelm 
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et al ., 199^ . 199 7) spectrum of the quiet solar corona, and are more consistent with a higher Alfven 
wave amplit ude, 11 0 km s. The remaining references give abundance measured in solar flares 


with RESIK (Sylwester et all . 120051 ). All models assume a constant amplitude with height for slow 
mode waves propagating up from the photosphere of 6.25 km s 1 . Fractionations are somewhat 
sensitive to this, in that increased slow wave wave amplitude decreases the degree of fractionation. 
However for the case considered here, that of a coronal Alfven wave on resonance with the coronal 
loop, most of the fractionation is restricted to a small region close to the top of the chromosphere 
w here the slow mode wa ves ge nerated parametrically by the Alfven wave driver are also significant. 
Rakowski and Laming] ( 2012 ) extended such considerations to Alfven waves away from the loop 
resonance. A similar FIP effect is produced, but the fractionation occurs over a wider layer of the 
chromosphere, as Alfven waves penetrate deeper. Some subtle differences occur. S is more strongly 
fractionated in such a case, and He is depleted much less, both as a result of the difference in Alfven 
wave behaviour. Slow mode waves are also more strongly generated deeper in th e chr omo sphere , 
and in cases such as this, their generation has a stronger effect on t he fractionatio n. 


Laming ( 2004a ) 


also considered the effect of different chro msopheric models f r om lYernazza et all (|198ll ). If VALC 
may be considered as a forerunner to the Avrett, andJLoese^ J 20081) model used here, and it does 


indeed give similar FIP fractionation, the result of iLaminpl d2004at) in that stronger fractionations 
are found within cells at dark points (VALA, VALB) and weaker fractionation at network segments 
(VALD, VALE). However these chromospheric models should probably be considered obsolete at 
the current time of writing. 


7.2 Open field 

We next consider the FIP effect in open magnetic field. The chromospheric model is extrapolated 
upwards assuming a density scale hei ght of 170 .000 km, to give a density profile comparable to 
that observed and modeled elsewhere (Laming], 2004b]). The temperature is set at a maximum of 
10 6 K, but is similarly inconsequential as in the closed field case. Figure [9] shows the same panels 
as Figure [TJ but extending to an altitude of 500000 km. Five Alfven waves are included, designed 
to match the spectrum given in Figure 3 of lCranmer et all ( 2007 ). with amplitude in the transition 


region from their Figure 9. Waves of angular frequency 0.010, 0.031, 0.062, 0.093, and 0.124 s“ 
are included with amplitudes at 500 000 km of 12.5, 150, 75, 50, and 12.5 km s -1 , respectively. The 
integ ration of Equations dill) begins at 500 000 km with outgoing waves only following Cranmer 


et al. ( 20071 ). and is taken back to the chromosphere. No account is taken of the motion of the 
coronal hole plasma. For the range of altitudes considered, the fast solar wind speed is always 
much less than the Alfven speed, and is probably still insignificant compared to the uncertainty in 
the Alfven speed. 

Figure [TUI shows a similar set of panels to Figure [8] above. Top left shows 5v and 
for the five waves in the chromosphere, with the real and imaginary parts distinguished as before. 
Bottom left now shows the chromospheric wave energy fluxes. Although we have made no assump¬ 
tion about where the waves originate, in the case of the coronal hole the usual interpretation is 
that Alfven waves start out as kink mode waves on photospheric flux tubes, and propagate up 
into the coronal hole. This bottom left panel shows that most of the upgoing wave energy flux is 
reflected back down again, at the steep chromospheric density gradient. Of order 10% of the wave 
energy makes it out of the chromosphere, and more reflection occurs higher up (not included in 
our model). Higher frequency waves suffer less reflection than the low frequency waves. 

The top right panel shows the ponderomotive acceleration. Superficially it is very similar to 
the previous case, with a “spike” at the altitude where the chromospheric density gradient is 
steepest. However, the strong wave reflection leads to a stronger component of the ponderomotive 
acceleration at lower altitudes, around 1000 km above the photosphere, and a reduced accelera¬ 
tion associated with the “spike”. In principle, stronger fractionation can occur low down in the 
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Figure 9: Solution for the coronal hole. The top panel shows <5v and <56/^47rp (Equations (fT2l) h 
real components in black, imaginary components in grey. The middle panel shows the wave energy 
fluxes, left and right going, and their difference, divided by the magnetic flux density. Normalized 
thus, this difference should be constant in the absence of wave growth or damping, and provide 
a check on energy conservation in the integration. The bottom panel shows the ponderomotive 
acceleration, a. 
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Figure 10: The top and bottom panels on the left hand side and the top panel on the right hand side 
show the same variables as in Figure [9] Additionally, the amplitude of slow mode wave generated 
parametrically by the Alfven waves is also shown in the top right panel. Compared to the closed 
loop case with a wave on resonance, the ponderomotive acceleration has a stronger contribution 
lower down in the chromosphere. The bottom right hand side panel shows the chromospheric 
ionization balance and FIP fractionation for selected elements. The bottom panel shows FIP 
fractionation for the abundance ratios He/O, Mg/O, and Fe/O, with the linestyles indicated as 
black lines, to be read on the left axis. The degree of fractionation is reduced from the closed field 
case. The ionization balance is shown as grey lines with the same linestyles for Fe, Mg, Si, S, and 
He, and is very similar to the previous, except that in the reduced radiation field of the coronal 
hole, He stays neutral longer. 
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chromosphere because here the hydrogen is mainly neutral, whereas around the “spike” it is be¬ 
coming ionized. This difference in the background H makes a big difference to the value of u e s in 
Equation m, and also to the resulting fractionation. The degree to which this actually happens 
is determined by the amplitude of slow mode waves appearing in v1 in the denominator of the 
integrand in Equation (12211 . The model shown in Figure [TO] assumes the same slow mode wave 
amplitude as in the closed loop case, a constant value of 6.25 km s _1 , which when added in quadra¬ 
ture with the parametrically generated slow mode wave amplitude eliminates the fractionation low 
down. Thus the fractionation is restricted to the top of the chromosphere, and generally agrees 
well with coronal hole observations. Treating t he Alfv en a nd acoustic waves in a more complete 
fashion (but only considering Fe/O). ICranmer et all ( 200 7) reach a similar result. The amplitude 
of acoustic waves developing low in the chromosphere in their model is higher than in ours, of order 
10 km s- 1 from their Figure 9, though varying with altitude. 


Table 4: FIP Fractionations in Open Magnetic Field (see text for details). 


ratio 

models 
6.15 6.25 

( km s -1 

6.5 

) 

a 

observations 

b 

c 

He/O 

0.68 

0.68 

0.68 

0.59-0.63 

0.55-0.69 

0.45-0.55 

C/O 

3.69 

1.00 

0.99 

1.41 -1.68 

1.41-1.68 

0.9-1.1 

N/O 

0.89 

0.89 

0.89 

1.07-1.32 

1.07-1.32 


Ne/O 

0.81 

0.81 

0.81 

0.44-0.52 

0.44-0.52 

0.3-0.4 

Na/O 

7.12 

1.56 

1.43 


1.45-1.91 


Mg/O 

6.74 

1.48 

1.37 

1.61 1.85 

1.29-2.82 

0.95-2.45 

Al/O 

6.46 

1.41 

1.33 


1.51-2.00 


Si/O 

5.97 

1.31 

1.26 

1.86-2.26 

1.29-2.34 

0.9- 1.8 

P/O 

5.22 

1.18 

1.15 




S/0 

4.80 

1.13 

1.11 

1.46-1.60 

1.17-1.86 


Cl/O 

1.01 

1.00 

1.00 




Ar/O 

0.95 

0.95 

0.95 




K/O 

7.76 

1.62 

1.48 




Ca/O 

7.77 

1.62 

1.48 


1.38-1.82 


Ti/O 

7.83 

1.62 

1.48 




Cr/O 

7.78 

1.60 

1.47 


1.81-2.63 


Fe/O 

7.80 

1.60 

1.47 

1.45 -1.80 

1.51-2.29 

0.65-1.35 

Ni/O 

7.25 

1.50 

1.39 




Kr/O 

0.98 

0.98 

0.98 





The effect of varying the slow mode wave amplitude is illustrated in Table [U where models are 
given for amplitudes of 6.15, 6.25 and 6.5 km s _1 . Above about 6.25 km s _1 , the FIP fractionation 
slowly decreases with increasing slow mode waves, assuming the Alfven waves are kept constant. 
Bel ow this value, the FIP fr actionation increases dra matically. Observa tional ratios are taken from, 
(a) IZurbuchen et all (2002h . given rel ative to O, (b) Bochslerl (2007a 1. relative to O. and (cl Ko 


et al. (120061 ), relative to H. FIP fracti onations at the level of those in the first mode l colum n in 
Table |H or higher, have been reported ( Widing and Feldman L 1992^^3mig_and_Mason . 1997ll. and 


Del Zanna 


are a l ways observed in open f ield structures. Other authors ( Doschek and Lamind . 2000 
120031 : iDel Zanna et all [20031 1 ha ve not found large FIP effect s in polar plumes. This last reference 


actua lly challenges the finding of Widing and Feldman ( 1992t h (but not that of Young and Mason 
1997th arguing that an isothermal plasma at logT = 5.9 would produce the same Ne VI/Mg VI 
intensity ratio with photospheric abundances. Hence the reality, and certainly the ubiquity of 
strong FIP effects in polar plumes is questionable. Even so, the possibility of strong FIP effects in 
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these structure is supported by our model, as Alfven waves can have considerable amplitudes low 
in the chromosphere, and the sensitivity of our open field FIP model to the assumed slow mode 
wave amplitude is perhaps realistic. 


7.3 The helium abundance 


The ponderomotive force was originally invoked to explain just the FIP effect ( Lamin d . 2004a 1. 
i.e., the enhancement in coronal abundance of the low FIP ions. Once more accurate treatment s 
involving the non-WKB analysis of Alfven wave transport were implemented (Laming, 20091 . 2012 1. 
it was noticed that as well as providing a satisfactory explanation of the FIP effect, the model also 
predicted a depletion of the He abundance, (and to a lesser extent Ne also) with respect to O in 
the corona and solar wind. As noted in Section [3j in the slow solar wind especially, the He/H ratio 
is reduced from its photospheric value and becomes more variable in slower speed slow solar wind. 
The He abundance is also depleted in the fast solar wind, but to a lesser extent, and is much less 
variable. 


Rakowski and Laming ( 20121 ) investigated models for the depletion of He in more detail. They 
showed that similar behavior is seen in the ratio He/O measured by ACE/SWICS (see Figure [IT] left 
panel); a depletion of He to about 0.4-0.8 of the solar photospheric value, with greater depletion 
seen the slower the solar wind speed, although th e sola r cycle dependence is not apparent in the 
time interval 1998-2011. Rakowski and Laming (l2012h constructed mode ls for a variety of loop 
lengths and corona magnetic fields, all using the Avrett and Loeseil (2008) model chromosphere, 
and studied the fractionation produced as a function of the assumed Alfven wave frequency with 
respect to the loop resonance frequency. For each case, the coronal wave amplitude was chosen to 
give a fractionation Fe/O by a factor of 4, as observed in the slow speed solar wind. As we have 
seen above, for waves that resonate with the coronal loop, most of this fractionation occurs at the 
top of the chromosphere, whereas off resonance, (or in open field regions) fractionation lower down 
in the chromosphere is possible. 



Figure 11: Left: He/O relative to the solar photospheric ratio as a function of the wave frequency, 
in units of the loop resonance frequency, for different loop lengths and magnetic fields. A trend of 
increasing He depletion with decreasing magnetic field strength, and possibly also increasing loop 
length, can be seen. Right: Wave energy transmission coefficients for the same loop models as a 
function of wave frequency. Shorter wavelengths (lower magnetic field, shorter loops) have broader 
resonances an^arejnore^^easil^^ransmittedAhe^^chromosjAiere^coronaJ^oundar^^ron^Rakowski 
and Laming (l2012ll . Reproduced bv permission of the A AS. 

Results are shown in Figure [TT] left panel. Helium is seen to be strongly depleted with respect 
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to oxygen when the wave frequency coincides with the loop resonance. In these cases, the pon- 
deromotive acceleration is restricted to the top of the chromosphere. This means that O can be 
accelerated into the corona once it becomes ionized, while He, the last element to remain neutral, is 
left behind. Elsewhere, where waves are not on resonance, the ponderomotive acceleration develops 
over a greater range of chromospheric altitudes. Correspondingly the acceleration in the region 
where O is ionized is smaller, since we are restricting ourselves to cases where Fe/O fractionates 
to a factor of 4. Hence O is not accelerated relative to He to the same extent, and the He/O ratio 
is unaffected, while otherwise the usual FIP effect still develops, albeit with subtle changes in the 
fractionation pattern. 

Closer inspection of the left panel of Figure [TT] also reveals greater depletion of He at the 
resonance with weaker magnetic fields. Again, the ponderomotive acceleration is becoming more 
concentrated towards the top of the chromosphere, but for a different reason. As the magnetic field 
weakens, the chromospheric layer where the plasma /? = 1, the equipartition layer, moves to higher 
altitudes. At this layer, a myriad of wave phenomena occur; reflection, transmission and mode 
conversion, and we do not attempt to model this except to say that FIP fractionation must occur 
above this layer. Hence as this layer moves upwards, ponderomotive acceleration is more restricted 
to the upper chromospheric layers, and He/O becomes more depleted. This is most easily seen in 
Figure QT] for frequencies just above the resonant frequency. 

In principle, the He/O depletion should also depend on loop length. Longer wavelength waves, 
resonating with longer loops, are more effectively reflected at density gradients. Thus we would 
expect more He/O depletion to be associated with longer coronal loops, and this seems to be borne 
out by the three loop lengths with magnetic fields of 15G. The three cases with 10G coronal fields 
are less clear, mainly because for many wave frequencies, and Fe/O enhancement of 4 could not 
be achieved, and so these points are not plotted. The increased penetration of shorter wavelength 
waves into a density gradient before reflection is illustrated in the right panel of Figure [Tl] Wave en¬ 
ergy transmission coefficients for the same loop models are plotted as a function of wave frequency. 
Waves with shorter wavelengths (lower magnetic field, shorter loops) have broader resonances and 
are more easily transmitted across the chromosphere-corona boundary. 

The fact that lower values of He/O are to be found for longer loops and weaker magnetic fields 
may have some bearing on the origin of the slow speed solar wi nd, and more specifically the origin of 
its different speed components. Some authors (e.g., Cranmer et all 2007]) argue that both fast and 
slow speed solar wind originate in open magnetic flux tubes, and that the difference in fractionation 
arise from subtle differences in the wave propagation on such structures. We would argue here that 
different slow mode wave amplitudes lower down in the chromosphere are the most likely variable. 
But in this case, the He depletion in fast and slow should be similar. Observationally this is clearly 
not the case, and our inference that the He/O depletion is strongly dependent on the Alfven wave 
frequency with respect to the loop resonance, being strongest when on resonance suggests that 
the slow solar wind plasma must originate in closed magnetic loops where it becomes fractionated, 
and that the Alfven wave causing the fractionation must be generated in the loops themselves, in 
order that the loop resonant frequency is selected. 


7.4 Significance of coronal Alfven waves 

The Alfven^aveJnvel£_su£gestedmbove_are_simila^o_those_observedAn_soln^ares_(e i g^Alexan- 


der et a,L. 1199811. but higher than those usually considered and observ ed in the solar corona (cf. 
Peteil . 200 lT De Pontieu et ai . 2007 : Peter, 201(ll : McIntosh et fl/ll201ll ). The predicted slow mode 
wave amplitudes actually match qu ite we ll with the non-thermal broadening observed close to loop 
footpoints in an active region bv [Baker et al. ( 2013 1. and correlated with the locations of FIP 
fractionation (see Fig. El. The Alfven wave values can be quite reasonable if the MHD fluctuation 
is confined to a small part of the coronal flux tube cross section. In this case a two component 


41 










































line profile should be expected, as in Peter (2001), with a narrow component with nonthermal line 
broadening with peak amplitude of ~ 25 km s _1 , and a broad component corresponding to Alfven 
waves with peak amplitude of ~ 100 km s _1 . Depending on the filling factor of the strongly oscil¬ 
latory plasma, this second component may or may not be readily detectable. Fila mentar y m odels 
of co r onal and flare heating have been invoked to explain flare lightcurves (Warren and Doschek, 
2005; IWarrenl . I200(tI 1 and stellar coronal emission measure distributions ( Cargill and Klimchuk , 


2006s). Figure [12] shows in the left panel the results of a 3D compressible MHD simulation includ¬ 


ing parallel heat conduction and radiation , of a coronal loop subject to forcing at its footpoints 
by photospheric motions (j Dahlbu rg et all 2012|). The plane depicted is across the midpoint of 


the loop, and one can easily see localized temperature hotspots developing. Physically, the mag¬ 
netic field is becoming stressed and periodically releasing the stored magnetic energy in small scale 
reconnection event. The right panel show observations with the Hi-C instrument (iTesta et al. 


20131 ) of the footpoint region of a coronal loop. Filamented, and variable emission in Fe XII can 


be see, interpreted as the thermal conductive response to filamented and sporadic heating events 
higher up in the coronal portion of the loop. The left panel is scaled to 4 000 km on a side, and 


so represents a similar spatial scale to that observed by Testa et al. (2013). The Fe XII emission 


peaks at a temperature of about 1.6 x 10 6 K, just above the top of the temperature scale in the 
simulation. The filamentation appears on an angular scale of arseconds or smaller, and thus is only 
resolvable by Hi-C. Scales of this size would have been below the angular resolution of prior imag¬ 
ing spectrometers or imagers, and consequently any Alfven waves associated with such structures 
would have been difficult to detect. 



0 2 4 6 8 10 12 0 2 4 6 8 10 12 


x [arcsec] x [arcsec] 


Figure 12: Left: Results of 3D MHD simulation from lDahlburg et al ( 2012 1. Temperature contours 
in the loop mid-plane are shown, illustrating the filamentary nature of the loop heating. The 
length scale is 4000 km on a side. The temperature scale runs from 4.9 x 10 5 K to 1.3 x 10 6 K. 


Reproduced by permission of the ESO. Right: Observations of active region moss from ITesta et ali 
( 2013t f. Filamentary and variable emission is seen, most likely as a response to filamentary heating 
in the coronal portion of the loop, communicated to lower altitudes by heat conduction. At 1 
arcsec ~ 750 km, this structure is of comparable size to the simulation. Reproduced by permission 
of the A AS. 


The filamentation is usually taken to derive from nanoflare heating (see below), though it 
could also refer to models of coronal heating by Alfven resonance. In a flux tube with cross-B 
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density gradient scale length l, a kink mode oscillation of the flux tube becomes Alfvenic at a 
resonant surface where the wave frequen c y u j = 2 L/Va, where L is the loop length and Va is 
the Alfven speed (see, e.g., iRuderman and Roberts! . l2002h . The width of this resonant layer is 
5 ~ (Iv/lo) ' , where v = oVa/R is the kinematic viscosity, a is the loop radius and R is the 
Reynolds number. The Alfvenic velocity fluctuations are larger than those of the kink mode by 
a factor 1/5 = ( l 2 R/aL ) 1 ^ 3 ~ O.li? 1 / 3 under typical conditions. Such wave motions observed in 
the solar corona have been identified as kink fa subset of fast model mode waves (e.g.. Nakariakov 
et aL l19991 : 1Wang and Solankil . l2004l h Alfven w aves w hich ge nerate the ponde romotive force in the 
chromosphere have been much harde r to detect rfV an Doorsselaere e t all 120081b though see Erdelvi 


and Fedun d2007t l and iTomczvk et all ( 2007t >. The global kink modes above presumably derive 
from acoustic motions in the chromosphere driving loop footpoints initiating the loop oscillation, 
that then subsequently decay into Alfven waves at the resonant surface. Turbulence associated 
with heating on particular field lines will only produce fractionation if chromospheric upflows are 
restricted to these field lines. However, if the upflow is the result of chromospheric evaporation, 
this is precisely what we should expect. 

The nanoflare paradigm suggests that the dominant loop footpoint motions are of much lower 
frequency, and do not excite oscillations in the loop but act so as to build up magnetic stresses 
in the corona. These stresses periodical ly release themselves, in what has become known as a 


“nanoflare”, as a current sheet develops ( Parker! . 1988 b Rannazzo et al\ 1 2007 . 12008 1 discuss the 
buildup of magnetic stresses within the framework of turbulence phenomenology, where it appears 
that veloci ty perturbations similar to the ~ 30 km s - 1 observed should be expected. Longcope 


et al. ( 2009h conjecture that in impulsive reconnection in post-flare loops, only about 10% of the 
liberated magnetic energy is converted directly into heat, the rest reappearing as kinetic energy that 
ultimately drives turbulence. In the case of reconnection of field lines at angle 6 1 equating magnetic 
energy destroyed to kinetic energy gained, B 2 sin 2 (9/2) /87r = pv 2 / 2, suggests v = va sin0/2. For 
B = 20 G and p = 1.67 x 10 -24 x 10 9 g cm -3 , v^_ = 1400 km s -1 , and y = 100 km s _1 implies 
0 = 8°, slightly larger than the expectation of lRappazzo and Parker ( 2013 b Of course some 
magnetic energy may go directly to heat, but it is plausible that such wave generation explains 
why surveys to find locali zed hot plasm a as evidence of nanoflare reconnection have generally been 
unsuccessful (e.g., Warren et all l20lH b Instead, energy goes from magnetic field to waves, and 
is thus gradually dissipated as heat throughout the corona, and not quickly and locally as might 
have been exp ected . 

Sturrock (|1999T ) gives a pedagogic review of the mechanisms by which various wave modes may 
be excited by reconnection. The reconnected field line is generally distorted, and this can either 
propagate away from the reconnection site as an Alfven wave, or em i t ma gnetoacoustic waves 
traveling perpendicularly to the magnetic field direction. Ilsobe et~. (20081) model small scale 
reconnection in the chromosphere. The emerging magnetic flux reconnects with the previously open 
held, to produce a reconnection jet accompanied by an upward propagating Alfven wave, which 
appears to be of appropriate frequency (0.01 Hz) and amplitude (20 km s _1 in the transition region) 
to give rise to some fractionation. Observational evidence of low-lying reconnecti on p ro ducing jets 
and transverse waves of similar frequencies and amplitu d es has been r e porte d ( Nishizuka et a,l\ . 
20081 : He et a,L . 20091 : Vasheghani Farahani et a,L . l2009h . Kigure et all ( 201fll ) explicitly consider 
the generation of Alfven waves by magnetic reconnection, and find that a significant fraction of 
the magnetic energy released (several tens of %, depending on geometry and plasma /3) can be 
carried off by Al fven or magnetoac oustic (fas t or slow mode) waves, with Alfven waves dominating 
for j3 < 1. Liu et al. (2 01 l aO and lLiu et al\ ( 201 lbl ) discuss the role of temperature anisotropies 
and wave_generation by the firehose instability in the outflow. 

Van Ballegooijen et al. ( 201 il l challenge some of these ideas and offer a more traditional view of 
an Alfven wave heated corona, where Alfven waves are introd uced at loop footpoints from below and 
are either transmitted or reflected, following iHollweg ( 1984l h They argue that “nanohare” heating, 
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as in Parker’s concep t, cannot supply sufficient heat to the corona. However, the sim ulations of 
Dahlburg et al. ( 2012h . (see also Gudiksen and Nordlundl . [2004lBingert and Peter . 2011 1. discussed 
above seem to contradict this statement. A million degree corona can indeed be heated and 
maintain ed by r ando m m otions of loop footpoints with timescale longer than the loop resonance. 
Further, Rakowski and La ming! (2012 ) show that ex tended Alfven wave propagation throughout 
the chromosphere, as in van Ballegooiien et SB, is unlikely to correctly predict the depletion 
in the coronal abundance of helium. 




Inclination (radians) 


Figure 13: Left: Schematic diagram from Khomenko and Callvi ( 2012 ) of acoustic to magnetosonic 
wave mode conversion at the plasma /3 = 1 layer of the chromosphere. An upcoming acoustic 
(fast mode) wave where va cs (/? ^ 1) can be transmitted into the overlying region where 
va -C cs (/? <C 1) as an acoustic (slow mode) wave, or mode converted to a magnetoacoustic 
(fast mode) wave. The magnetoacoustic wave refracts in the region where va is increasing with 
height and eventually rejoins the lower (3 1 plasma. The transmitted slow mode can continue 

propagating until it reaches an acoustic cut-off. When the plane in which the fast mode refracts 
is distinct from the plane of inclination of the magnetic field (indicated here by blue lines), it may 
also mode convert to an Alfven wave, which is not considered further in this paper. Reproduced 
by permission of the A AS. Right: Transmission coefficients for isotropic acoustic waves (in the 
upgoing hemisphere) as a function of the inclination of the magnetic field, for different values of 
kHr>, calculating by integrating a generalization of Equation (l30l) in the region O<0<7r//2 and 
0 < (f> < 27t. 


7.5 Modeling the inverse FIP effect 

One of the striking features of the Wood-Linsky relation illustrated in FigureS] is how the coronal 
abundance anomaly smoothly changes from solar-like FIP bias at spectral types G to early K, to 
inverse FIP in M dwarfs. This suggests that a model for solar-like FIP fractionation should, with 
suitably chosen parameters, be capable of predic t ing an inverse FIP effect. 

The model described above and in lLamin 5 ( 2004^ . 20091 . 12012 4 for the FIP effect suggests 
that Alfven waves of amplitude approximately 100 km s _1 are generated in the coronal loop at 
the resonant frequency, and remain trapped in the loop “resonant cavity”. Upon reflection from 
chromospheric footpoints, the waves develop a ponderomotive force in the steep density gradients 
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there, and this force, acting on chrom osp heric ions (but not neutrals), preferentially accelerates 
these ions up into the coro na. Laming ( 2012) st udies the fractionations produced by waves on and 
off resonance, and iRakowski and LamingT i 2012 1 extend this to different loop lengths and magnetic 
fields, concentrating mainly on the fractionation of He with respect to O. Th ese works onl y consider 
cor onal Alfven waves, with chromospheric acoustic waves included following iHeggland e.t al. ( 2011 1 
and lCranmer et aZl ( 2007h as terms in the denominator of the integrand in Equation (|22|) . When 
these upcoming chromospheric acoustic waves are allowed to mode convert, at the layer where 
sound and Alfven speeds are equal, to what in the magnetically dominated upper chromospheric 
become fast mode waves, inverse FIP fractionation can result. This arises because the fast mode 
waves undergo reflection back downwards as the Alfven speed increases, giving rise to a downwards 
directed ponderomotive force than can compete with that due to the coronal Alfven waves. Inverse 
FIP requires d5E 2 /dz <0. If waves propagate from chromosphere to corona (or from corona to 
chromosphere) without reflection, then dSE 2 /dz > 0 always. This is because SE = SvB (where 
B is ambient field, assumed uniform here) and 8v increases as the density decreases. If waves 
propagating up from beneath are reflected back down again, Sv and SE must decrease with height, 
and dSE 2 /dz < 0, giving rise to a downward ponderomotive force. 

The p r ocess of m o de con version is shown schematically in the left panel of Figure ll3l taken from 
Khomenko and Callv ( 2012l l. An upcoming acoustic (fast mode) wave where va "C cs (/3 1) can 

be transmitted into the overlying region where va cs {/3 -C 1) as an acoustic (slow mode) wave, 
or mode converted to a magnetoacoustic (fast mode) wave. The magnetoacoustic wave refracts 
in the region where va is increasing with height and eventually rejoins the lower /3 1 plasma. 

The transmitted slow mode can continue propagating until it reaches an acoustic cut-off. When 
the plane in which the fast mode refracts is distinct from the plane of inclination of the magnetic 
field (indicated here by blue lines), it may also further mode convert to an Alfven wave, which is 
not considered f urt her in_this paper. The acoustic waves have an energy transmission coefficient 
of (Callv and Goossens . 2008 1 


T = 


exp 


—7T |k| sec 0^1 — (sin a cos </> sin 9 + cos a cos 9 ) 2 ^j 


[d(V 2 /c? s )/dz ] p=1 

= exp 7T |k| Hd seed ^1 — (sin a cos <j> sin 9 + cos a cos 0) 2 ^ 
—> exp (—7r |k| Hd sec 9 sin 2 9) , 


(30) 


where a and 9 are the angles to the vertical made by the magnetic field and the wavevector 
respectively, and (j) is the polar angle between them. The last line gives the more familiar result 
in vertical magnetic field. The right panel of Figure [13] shows the acoustic to acoustic energy 
transmission coefficient for initially isotropic waves (in the upward going hemisphere) as a function 
of magnetic field inclination to the vertical (the angle a above), for values of |k| Hr) = 0.25, 0.5, 
and 1. For solar oscillations with period 300 seconds, kHj) = ^/w 2 Ff|,/c| — 1/4 = 0.25 where 
we have taken Hd = 168 km from Figure [6] and cs = kT/m p cs 6.3 km s . In later type stars 
with deeper convection zones, we might expect kHo to increase, si nce u> ~ q/cs incre ases with 
increasing g and d ecreasing ec in the colder photospheres (see. e.g.JBnmtt et all 120101: Kicldsen 


and Bedding, 201 1), while Hd in the low chromosphere remains approximately constant, as far as 
can be determined from model chromospheres ava i lable in the literature fe.g.. iFuhrmeister et all 
20051 : iHoudebine and Stempelsl . 1997 : Vievtes et al. I . 2005l h These factors would increase the degree 
of mode conversion in later type stars. Solar p-mode s are also known to decrease in intens i ty an d 
increase in width from solar minimum to maximum ( Chaplin et al. , 2000l : Simoniello et oil l2010h . 
This suggests that the trapping of acoustic modes within th e solar env elope becomes less effective 
at solar maximum, and allows energy to l eak out ([Pinter et al. . 200ll ). Sunspots have long been 
known to be sinks of p-mode energy (see Braun . 19951 . and references therein), most likely through 
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mode conversion or resonant absorption to Alfven or fast m ode waves. P-nr odes ar e also known 
to decrease in intensity for later spectral types, modeled by Kieldsen and Bedding (2011) as due 
to extra “leakage”. 

A further reason for the transition from FIP to inverse could be that in stars on the left hand 
side of Figure IH upcoming acoustic waves encounter an acoustic cut-off before reaching the /3 = 1 
layer, and never have a chance to mode convert. In the stronger magnetic fields of later type 
stars, the /3 = 1 layer is deeper in the chromosphere, and upcoming acoustic waves will reach this 
first, before any acoustic cut off, and mode convert to fast mode waves. The fast mode waves are 
immune to the acoustic cut off, but ultimately refract back downwards. 

We tr eat the fast mode waves as approximately isotropic in the upward moving hemisphere, 
following IWood et al. ( 2012 ). Then the fraction reflected at chromospheric height z is 


f R (z) — a /1 


' (A/3=l 


V 2 (z) + c 2 (z) 


(31) 


where zp—\ is the chromospheric height where mode conversion occurs. The ponderomotive accel¬ 
eration due to fast modes waves is then 


c 2 d 
4 dz 


SE 2 


5v 2 1 dSv 


Sv 2 df R 
4 dz 


(32) 


The two terms represent an upwards contribution arising as the fast mode waves increase in 
amplitude as they propagate through lower density plasma, and the downwards contribution arising 
from fast mode wave reflection. Evaluating 


df R 

dz 


c| (zp—i) V A dV A _ c| (zp= i) Vj 


(V 2 +c 2 ) 2 f R dz 


and assuming from the WKB approximation 
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1 d5v 
Sv dz 2 Hr 


4H d ’ 


(33) 


(34) 


where Hr, and Hr are the signed density and magnetic field scale heights, we find from equation 
5 


a = 


(Jr - l) - 
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l 


8 H d AH b 


r 2 

C S 


(*/9=l) 


(VI 


+ 4 ) 


8 H d \H b 



(35) 


Remembering that both Hr, and Hr are negative, and f R < 1, the first term in curly brackets 
is negative, giving rise to Inverse FIP effect, and the second term is positive, giving the more 
usual FIP effect. In conditions where Va cs, an overall downwards pointed ponderomotive 
acceleration requires \Hr>\ < \Hb \ /6 in this simple model. 

Additional reflection of fast mode wave from, e.g., density fluctuations (not included in this 
model) would relax the requirement. So too would fast mode waves spreading out laterally from a 
horizontally localized source. Even so, Equation [35] implies that inverse FIP effect is more likely to 
found in stars with minimal magnetic field expansion through the chromosphere, which fits with 
its observation in M dwarfs. While the magn etic fields measured in these stars are similar to thos e 
in the Sun, the filling factor is higher (e.g., Donati and Landstreet . 20091 : Reiners e.t al . 2009l f. 
allowing less volume for expansion with increasing altitude. 

Figure mi illustrates a cal culation desi gned to give an inverse FIP effect. The model is similar 
to that shown in Figure 3 of Lamind ( 2012 1. A loop of length 100 000 km, with a 80 G magnetic 
field is considered, with a resonant Alfven wave of frequency 0.28 rad s -1 . The magnetic field is 
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Figure 14: Illustration of a model demonstrating the origin of the Inverse FIP Effect. The top left 
panel shows the variables Sv and SB/y/inp for the coronal Alfven wave. The bottom left panel 
shows the upgoing (dashed curve) and downgoing (solid curve) wave energy fluxes. The dotted line 
shows their difference. The top right panel shows the ponderomotive acceleration. The positive 
contribution in the upper chromosphere (solid curve) comes from the coronal Alfven waves. The 
negative contribution lower down (dashed curve) comes from the total internal reflection of fast 
mode waves. The positive fractionation induced by the Alfven waves is suppressed by the choice 
of evaporative flow speed, which is higher in the lower density regions of the chromosphere due to 
continuity. Suppression of fractionation low in the chromosphere by slow mode waves is reduced 
in this example, because the slow modes mode convert to fast modes. 

The dotted curve gives the amplitude of acoustic waves through the chromosphere. The bottom 
right panel gives the FIP fractionations for the ratios S/O, He/O, Mg/O and Fe/O. He/O remains 
unchan ged, but all ot hers d isplay an inverse FIP effect. The choice of an “active region” spectrum 
from Vernazza and Reeves (119781 1 has increased the ionization fraction of He, compared to the 
earlier examples. 
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compressed by a factor 0.5 through the chromosphere (previously it was 0.2), and a fast mode wave 
amplitude of 10 km s _1 is included at the /? = 1 layer, which is allo wed to propagat e and refract as 
described above. The model chromosphere derives from lAvrett and Loesed ( 2008ll . though future 
work should implement a model stellar chromosphere. The top left panel shows the variation of 
the perturbations Sv and SB/y/lnp associated with the coronal Alfven wave, and the bottom right 
panel shows the Alfven wave energy fluxes, both upward and downwards directed. The dotted 
line in the bottom right panel shows the difference in wave energy fluxes. The top right panel 
shows the ponderomotive acceleration. The positive contribution in the upper chromosphere (solid 
curve) comes from the coronal Alfven waves. The negative contribution lower down (dashed curve) 
comes from the total internal reflection of fast mode waves. The dotted curve gives the amplitude 
of acoustic waves through the chromosphere, modeled as outlined above. The bottom right panel 
gives the FIP fractionations for the ratios S/O, He/O, Mg/O and Fe/O. He/O remains unchanged, 
but all others display an inverse FIP effect, coming from the fast mode waves. The positive FIP 
that would arise from the Alfven waves higher in the chromosphere is suppressed by the upward 
flow speed through the chromosphere, taken here to be 10 6 cm s -1 at a chromospheric density of 
lO 10 cm -3 , and of course smaller lower down in the chromosphere due to the increased density. 
This upflow is now included in the Alfven wave transport equ ation 11. Up flows of at least this 
speed are typical of the gradual phase of solar flares (e.g., Czavkowska et ail 19991 ). High in 
the chromosphere, the He ionization fraction is i ncreased r elative to the earlier examples by the 
choice of an “active region” spectrum fron lVernazza and Reevesl (1978) with which to illuminate 
the chromosphere. This will reduce any depletion of He that might otherwise occur in conditions 
giving rise to positive FIP effect, and might be relevant to the relatively high He abund ances 
observed in flares and CMEs (e.g. Feldman et all 20051 : Wimmer-Schweingruber et all 20061 ). 

The transition from FIP effect to inverse FIP effect with increasing fast mode wave amplitude 
at the /? = 1 layer is illustrated in Table [5] Models with fast mo dewave amp litudes of 5, 10, and 
15 km s _1 are compared with abundances in M dwarfs taken from Liefke et al\ (20 081). who co mpare 
observed coronal abundances with solar photospheric a bundances of lAsnlund et al\ ( 2005d ). The 
coronal abundances for EV Lac given bv iLaming (2009|) are in better agreement with the models, 
especially for Ne/O and S/O, but even so, it is clear that the broad systematics of the inverse FIP 
effect are well reproduced by the mo dels in Table [5] Re sults for further minor ions observed in 
g Gem and HR 1099 are taken from iHuenemoerder et al. (2013a|b ). The Alfven wav e amplitude 


has been reduced from that in Figure [7] by a factor of two. Wood and Laming (2013) argue that 
this m ight be appropriate in the strong magnetic fields of late type stellar coronae: Drake et al .I 


( 20061 ) studied the efficiency of electron acceleration in reconnection with the ambient plasma f3 
(the electron plasma (3 e = 8Trn e kBT e /B 2 , to be more precise). They found maximum energy input 
to electrons at f3 e = 0, with reduced electron heating at higher /l e , or lower magnetic field. We 
suggest that at the left hand side of the Wood-Linsky relation, coronal reconnection primarily 
generates Alfven waves that end up causing positive FIP fractionation when they reflect from the 
chromosphere. As one moves to the right, to later spectral type, coronal reconnection puts more 
energy into electrons, and less into waves. Consequently, the positive FIP effect diminishes. 

While inevitably a highly simplified model of chromospheric wave processes, we take the po¬ 
tential for the ponderomotive force to explain the Inverse FIP effect at later spectral types, as 
well as the more usual FIP effect in the Sun, as a significant point in its favor. The Inverse FIP 
effect is also observed in many other more active stars, and we expect a similar explanation to hold 
there. However, the wave origin will be more complicated when additional processes due to stellar 
rotation, tidal interaction (e.g., RS CYn binaries) and accretion (e.g., T Tauri stars) are included. 
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Table 5: FIP and Inverse FIP Effects 


Ion 

Ampl. (km 
5 7.5 

s’ 1 ) 

10 

EQ Peg A 

from Liefke et al. 12008! 
YY Gem AU Mic EV Lac 

AD Leo 

from Lamine (2009) from Huenemoerder et al. 12013a! 
Prox. Cen EV lac <r Gem HR 1099 

N/O 

0.99 

0.99 

0.99 

1.00 

1.75 

2.08 

1.52 

1.67 

2.94 


2.88 

1.58 

Ne/O 

0.98 

0.98 

0.98 

2.7 

2.1 

2.31 

1.71 

2.33 

2.06 

0.96 

2.45 

2.45 

Na/O 

1.53 

0.87 

0.39 








0.78 

1.07 

Mg/O 

1.49 

0.83 

0.36 

0.52 

0.34 

0.28 

0.38 

0.49 

0.72 

0.23 

0.38 

0.32 

AI/O 

1.44 

0.75 

0.30 








0.63 

0.54 

Si/O 

1.37 

0.71 

0.28 

1.21 

0.87 

0.51 

0.85 

1.20 

1.80 

0.52 

0.44 

0.36 

S/O 

1.14 

0.57 

0.22 

1.52 

1.27 





0.86 

0.46 

0.41 

Ar/O 

0.98 

0.98 

0.98 








1.95 

1.48 

K/O 

1.57 

0.64 

0.18 








6.46 

< 4.07 

Ca/O 

1.57 

0.63 

0.17 








1.20 

0.93 

Fe/O 

1.55 

0.47 

0.09 

0.45 

0.29 

0.26 

0.33 

0.52 

0.46 

0.24 

0.38 

0.32 










7.6 Saturation 


We conclude this section with a simple discussion of what might limit the absolute magnitude of 
the FIP fractionation. Alfven waves of arbitrarily high amplitude will eventually erase the density 
gradients that give rise to the ponderomotive force. We estimate the ponderomotive acceleration 
at which this will occur, and hence a limit on how high this acceleration can go, as follows. 

The density in a gravitationally stratified atmosphere with a ponderomotive acceleration is 


P = ^ ex P (^~9 Z / V s + J & 1 'eff/i's/v*dzj 


(36) 


where all symbols are defined above in subsection [631 and here the element s of interest is H. Then 

v eff 
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We substitute a = —5 v 2 /8Hjj to find 
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(38) 


where Hjjq is the density scale length when a = 0, i.e. in the absence of the ponderomotive 
acceleration. Consequently 


a 0 


l _|_ £ v ‘ff fa 2 


1 


]__Hm _ 

a 0 v 2 1 + (1/^ - 1) v S i/v s 


(39) 


where ao = —8v 2 /8Hdq, the ponderomotive acceleration in the unmodified density gradient. Hence 
as ao —> oo, a —> — (1 + (1/^ — 1) v S ilv) v 2 /Hdq which is in general comparable to the ponderomo¬ 
tive accelerations invoked in this paper to explain the FIP effect, if we take £ = 1, with v s ~ 10 6 
cm s -1 and Hdo ~ 10 6 cm. A small departure from full ionization (^ < 1 increases this estimate 
quite quickly, because u S i » u sn . In this case we would have to appeal to properties of the Alfven 
waves themselves to explain the relative constancy of the FIP effect. If the waves are restricted to 
filamentary sections of a coronal loop, and are not monolithic oscillations of the whole loop itself, 
then regions of velocity shear will exist between oscillating and non-oscillating regions of the loop. 
This velocity may be expected to excite further wave motions, for example drift waves, which could 
provide another means of saturating the ponderomotive acceleration. Further developments along 
these lines are beyond the scope of this review. 


50 











8 Conclusions &; Future Work 


In this review we have attempted to show how the interpretation of abundance enhancements or 
depletions arising due to the action of the ponderomotive force due to various MHD waves offers the 
potential of a rather complete description of the phenomenon. It may also be hoped that advances 
along these lines, together with the detection and observation of waves in the solar atmosphere 
will be a profitable route towards solving the problem of coronal heating. Here we highlight some 
of the weak spots in the model, and areas where extra effort might provide a significant advance. 

1. Coronal Waves: Concentrating first on the solar FIP effect, our models above give a good 

account of the observed abundance anomaly with a coronal Alfven wave amplitude of around 
100 km s _1 . This is a peak wave amplitude, and should therefore be expected to give rise to 
a nonthermal line broadening around 70 km s _1 , measured at half maximum line intensity. 
This is higher than usually observed, but is not necessarily a problem if the line width is 
“diluted” by the filamentary structure across the loop. One important question is then to 
understand quantitatively where this velocity amplitude originates. We have sketched out 
some possibilities above, nanoflare associate d rec on necti on and Alfven resonance, but this 
is an area where numerical simulation (e.g., Dalilburg et all l2012h should be expected to 
yield new insights. The amplitude of waves generated will probably depend on the values 
of the coronal resistivity and viscosity, and thus have implications for theories of coronal 
heating. Different mechanisms of coronal heating might also be expected to produce Alfven 
waves of different polarizations. The examples given in this review have treated shear Alfven 
waves, such as might be expected to derive from reconnection. Alfven waves resulting from 
resonant absorption would be torsional waves, and those generated by a reconnection outflow 
by a firehose instability would exhibit some circular polarization. The Alfven wave transport 
equations and fractionation (Sections 16.41 and 16.511 do not change with polarization, but 
the coupling to other wave modes (Section 16.611 does. Pure parallel propagating circularly 
polarized waves do n ot cou pl e to slow m odes at all, while the case of torsional waves has 
been considered bv IVasheghani Farahani et al\ (|201 1 1. On twisted flux tubes, pure torsional 
Alfven waves do not exist, and some degree of mixing with the compressional kink mode is 
inevitable. It remains to be seen whether these changes in the wave physics will result in 
detectable changes in FIP fractionation, and if so, whether this represents a new avenue of 
approach to the problem of coronal heating. 

2. Chromospheric Waves: The solar FIP effect is subtly different depending on the altitude 

in the chromosphere where the fractionation occurs. High in the chromosphere, where H is 
becoming ionized, Fe and Mg fractionate more than Si, and S behaves more like a high FIP 
element. He can also be depleted relative to O. In the lower chromosphere where H is neutral, 
the low FIPs fractionate to essentially the same degree, and S behaves more like a low FIP 
element. He/O remains unchanged. Such fractionation can occur with off-resonant Alfven 
waves, e.g., in an open field region, or with upcoming fast mode waves in sufficiently diverging 
magnetic field. An important goal for spectroscopic and in-situ observations should be to 
obtain data of sufficient quality to distinguish between these tw o possi bili ties . Possibly the 
best extant spectroscopic analysis, that of quiet solar corona bv lBrvans et al\ ( 20091 ) (given 
as column “d” in Table [3]) strongly favors fractionation at the top of the chromosphere. 

Measurements of the S abundance made in situ in the slow speed solar wind often sh ow 


signifi cantly higher values th an those obtained from spectroscopy (e.g. iGiammanco et al. 


2007: Reisenfeld et al., 2007). The comparison between these two forms of measurement 


is very cl ear in Figure 1. of Sc hmelz et all ( 20121 ). Another element predicted to behave 
similarly ( Rakowski and Larninj . l2012t ) is C. Unfortunately, spectroscopic measurements of 
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the C abundance in the solar corona are very difficult, because at coronal temperatures, C 
is typically fully ionized and emits no lines. 


Fractionation low down in the chromosphere can occur with Alfven waves well away from 
the loop resonant frequency, such that they are not trapped in the coronal loop, or with 
upward propagating fast mode waves in region where the magnetic field expands, giving rise 
to positive FIP fracti onation from Equation (13511. He is w ell known to be present with high 
abundance in CMEs ( Wimmer-Schweinsrruber et al. , 20061 1 , consistent with fractionation low 
in the chromosphere, but S, the other element crucial in this regard is typically not measured. 
Whether fractionation occurs low down or high up depends not only on the coronal Alfven 
wave, but also on wave physics in the chromosphere. The mode conversion acoustic waves to 
fast mode waves at the /3 = 1 layer makes a connection with the fields of helioseismology and 
asteroseismology. There is much observational and theoretical work ongoing on these aspect 
of chromospheric wave propagation that is relevant to the FIP effect. One key component 
to this will be chromospheric vector magnetic fields. Currently chromospheric and coronal 
magnetic fields are extrapolated from photosplieric magnetic field observations. Uncertain 
knowledge of the magnetic field expansion through the chromosphere makes it difficult to 
know where mode conversion occurs, and whether upcoming fast mode waves should give 
FIP or inverse FIP effects. 


Considerable interest has arisen recently in the idea that the solar corona and wind might 
be s nnnlied bv Tvne II snicnles accelerated directly from the chromosphere (e.g.. De Pontieu 


et al ., 20 111 ; Martmez-Svkora et all 201 ill . rather than plasma being evaporated into a coro¬ 
nal loop before being released into the solar wind. We comment here that in such a case, 
it is difficult to see how a FIP effect could arise. With reference to Equation (1221) . the flow 
velocity associated with the spicule motion is likely to increase v1 in the denominator of the 
integrand to an extent such t hat the integral tends t o zero , and p s (z u ) = p s ( zi ), yielding no 
fractionation. In simulations, Martmez-Svkora et al\ 1 2011 1 find Type II spicule material ex¬ 
tending to tem peratures c l ose to 10 6 K, but no higher, consistent with the temperature range 
over which Laming et al. ( 19951 ) observed photospheric abundances (see Fig. [l]). Disconti¬ 
nuities in several other obs e rva bles have been located between 500,000 K and 10 6 K in the 
solar atmosphere ( Feldmanl . 1983, 1987), an d further consideratio ns of element abundances 
( Feldmanl . 1998ah reinforce the conclusion of Laming et al. (1995). Other authors have also 
questioned the conne ction between coronal mass supply and Type II spicules on t he basis of 


spectr al line profiles dPatsourakos et all 12014 iKlimchuk and Bradshawl . 120141 ) . [Goodman 


(120141 ) argues on energetics grounds that Type II spicules may power the quiet solar corona 
and coronal holes, but not active regions. 

Type II spicules are most likely the cause of the absence of the FIP effect at temperatures 
below about 10 6 K. However in the case of Procyon, with its lower surface gravity, Type II 
spicules are possibly the most likely reason why its corona appears to have photospheric 
composition, unlik e the o ther F dwarfs of similar spectral type in Figure [I] It is interesting 
to note that Drake et al. ( 1995al ). without knowing the cause, recognized that the FIP effect 
appeared in the full disk solar spectrum at about the temperature where the supergranula¬ 
tion disappeared from solar images, and speculated that if a similar transition happened on 
Procyon, it would occur at higher temperature. 


3. Solar Observations: Detailed spectroscopy and in situ measurements should aim to distin¬ 
guish between the two regimes of FIP fractionation mentioned above, i.e., low or high in 
the chromosphere. The aim should be to try and measure several element abundances si¬ 
multaneously, not just the evaluation of one or two abundance ratios. This will ultimately 
require spectrometers with bandpasses specially designed for the purposes. The Hinode/EIS 
instrument, for example, was designed in the 1990s (well before the ponderomotive force 
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model of the FIP effect was published) to maximize the coverage of lines from Fe ions, and 
so has rather few lines from high FIP ions within its bandpass, S being the most promi¬ 
nent. However, Hin ode/EIS does off er imaging spectroscopy sensitive to the FIP effect, as 
demonstrated by iBaker et all 1 20131 ). The observation of FIP enhanced regions correlated 
with high non-thermal broadening at loop footpoints is an important advance. The filamen¬ 
tary coronal heating, and correspondingly filamentary FIP fractionation could be potentially 
observable as a variation in FIP effect across the cross section of a coronal loop. Increased 
efforts to measure absolute abundances (i.e. relative to H) would also yield insights into other 
mechanisms working to modify element abundances, and be advantageous in comparisons of 
abundances from remotely sensed spectroscopic observations with those measured in situ in 
the solar wind. 


4. Stellar Observations: The quality of stellar observations is possibly less likely to improve in 
the future than those of solar observations, since missions in X-ray astronomy are fewer and 
further spaced in time than solar missions. Imaging spectroscopy is of course not possible, 
so the same comments about detailed spectroscopy to measure as many element abundances 
simultaneously apply. The main advance here will come from measuring abundances in a 
wider range of stellar targets. Observing more F dwarfs would flesh out the bottom left 
hand corner of Figure [4] and really establish a saturation of the FIP effect at about the 
level observed in the Sun, and no higher. This would then be interpreted as the FIP effect 
established solely by coronal Alfven waves, with no contribution from upcoming fast modes. 
Further observation with Chandra of r Bootis A would establish whether contamination from 
its M dwarf companion is responsible for its anomalous inverse FIP effect, or whether this is 
intrinsic to the F dwarf, and presumably due to the close-in Jupiter mass planet. 


Observations of inverse FIP effect exist in many more active and exotic stellar targets than 
those plotted in Figure [fl We have somewhat neglected discussion of these, for the simple 
reason that in trying to reach theoretical understanding of where the Inverse FIP effect comes 
from, it is obviously advantageous to start with simple objects first, before moving onto more 
complicated ones. Strong Inverse FIP seen in RSCVn binaries with tidal interaction would 
suggest that the planet hosted by r Boo A might produce the same effect. However, e 
Eridani also hosts a Jupiter mass planet, apparently with no unusual coronal abundance 
effects. While on the subject of stellar “complexity”, we note the poor state of knowledge 
of M dwarf chromospheres and photospheric abundances. This is an intrinsic issue, due to 
the complicated stellar spectra with many molecular bands and otherwise unidentified lines. 
There appears to be no “magic bullet” in sight, other than detailed careful work on obtaining 
and modeling such data. Further knowledge of aster oseism ology and associated wave physics 


may be expected to come from the Kepler mission ( Koch et al . 201 0l). 
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